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Abstract 

A derivation is given from first principles of the fact that the SU(2) gauge theory is in a 
confining phase for all values of the coupling < g < oo defined at lattice spacing (UV 
regulator) o, and space-time dimension d < 4. The strategy is to employ approximate RG 
decimation transformations of the potential moving type which give both upper and lower 
bounds on the partition function at each successive decimation step. By interpolation between 
these bounds an exact representation of the partition function is obtained on progressively 
coarser lattices. In the same manner, one obtains a representation of the partition function 
in the presence of external center flux. Under successive decimations the flow of the effective 
action in these representations is constrained by that in the upper and lower bounds which are 
easily explicitly computable. Conflning behavior for the vortex free energy order parameter 
(ratio of partition functions with and without external flux), hence 'area law' for the Wilson 
loop, is the result for any initial coupling. Keeping the string tension fixed determines the 
dependence g{a), which is such that g{a) — > for a ^ 0. 
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1 Introduction 



Four-dimensional SU{N) gauge theory at zero temperature is known to be in a confining 
phase for ah values of the bare coupling. A very large amount of work has been performed 
over the last decade in an effort to isolate the types of configurations in the functional 
measure responsible for maintaining one confining phase for arbitrarily weak coupling [I], [2]. 
Nevertheless, a direct derivation of this unique feature of SU{N) theories (shared only by 
non-abelian ferromagnetic spin systems in 2 dimensions) has remained elusive. 

The origin of the difficulty is clear. It is the multi-scale nature of the problem: passage 
from a short distance ordered regime, where weak coupling perturbation theory is applicable, 
to a long distance strongly coupled disordered regime, where confinement and other collec- 
tive phenomena emerge. Systems involving such dramatic change in physical behavior over 
different scales are hard to treat. Hydrodynamic turbulence, involving passage from lami- 
nar to turbulent flow, is another well-known example, which, in fact, shares some striking 
qualitative features with the conflning QCD vacuum. 

The natural framework for addressing the problem from first principles is a Wilsonian 
renormalization group (RG) block-spinning procedure bridging short to long scales. The use 
of lattice regularization, i.e. the framework of lattice gauge theory (LGT) [3], is virtually 
mandatory in this context. There is no other known usable non-perturbative formulation of 
gauge theory that gives the path integral in closed form preserving non-perturbative gauge 
invariance and positivity of the transfer matrix (reflection positivity). Attempts at exact 
blocking constructions towards the 'perfect action' along the Wilsonian renormalized trajec- 
tory [4], however, turn out, not surprisingly, to be exceedingly complicated. 

There are, nonetheless, approximate RG decimation procedures that can provide bounds 
on judicially chosen quantities. The basic idea in this paper is to obtain both upper and 
lower bounds for the partition function and the partition function in the presence of external 
center flux. The bounds are obtained by employing approximate decimations of the 'potential 
moving' type [5], [6], which can be explicitly computed to any accuracy by simple algebraic 
operations. This leads to a rather simple construction constraining the behavior of the exact 
partition functions in the presence and in the absence of center flux; and, through them, 
the exact vortex free energy order parameter. The latter is the ratio of these two partition 
functions. It is thus shown to exhibit conflning behavior for all values < P < oo, of the 
inverse coupling P = A/g'^ deflned at lattice spacing a (UV cutoff). An earlier outline of the 
argument was given in [7j. 

As it will become clear in the following, there are two main ingredients here that allow 
this type of result to be obtained. The flrst is the use of approximate decimations that 
are easily explicitly computable at every step, while correctly reflecting the nature of RG 
flow in the exact theory. The second is to consider only partition functions, or (differences 
of) free energies, rather than the RG evolution of a full effective action that would allow 
computation of any observable at different scales. This more narrowly focused approach 
results into tremendous simplification compared to a general RG blocking construction. 

The presentation is for the most part quite explicit. Some simple propositions, mostly 
containing basic bounds, serve as building blocks of the argument. They are enumerated by 
roman numerals in the text below. Most proofs have been relegated to a series of appendices 
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so as not to clutter what is essentially a simple construction. Only the case of gauge group 
SU{2) is considered explicitly here. The same development, however, can be applied to 
other groups, and, most particularly, to SU (3) which exhibits identical behavior under the 
approximate decimations. 

It will be helpful at this point to provide an outline of the steps in the argument devel- 
oped in the rest of the paper. In section [21 starting with the pure SU{2) LGT with partition 
function defined on a lattice of spacing a, we define a class of approximate decimation trans- 
formations to a coarser lattice of spacing ba. In section [3] the resulting partition function on 
this decimated lattice is shown to be an upper bound on the partition function on the original 
lattice. A similar rule can be devised for obtaining a partition function on the decimated 
lattice which gives a lower bound on the original partition function. One then interpolates 
between these bounds. For some appropriate value of the interpolating parameter, one thus 
obtains an exact integral representation of the original partition function. This representation 
is in terms of an effective action defined on the decimated lattice of spacing ab plus a bulk 
free energy contribution resulting from the blocking a — > ab. Now, any such interpolation is 
not unique, and it is indeed expedient to consider different interpolation parametrizations. 
The resulting partition function representation is then invariant under such parametrization 
variations in its effective action. The other important ingredient is that the effective action 
in this representation is constrained between the effective actions corresponding to the upper 
and lower bound partition functions. Iterating this procedure in successive decimations, a 
representation of the partition function is obtained on progressively coarser lattices of spacing 
a ^ ab ^ ab"^ ■ ■ ■ ^ aV^. 

In section [4] we consider the partition function in the presence of external center flux. This 
is the fiux of a center vortex, introduced by a Z{2) 'twist' in the action, and rendered topo- 
logically stable by winding around the lattice torus. The decimation-interpolation procedure 
just outlined for the partition function can be applied also in the presence of the external 
flux. A representation of the twisted partition function on progressively coarser lattices can 
then be obtained in a completely analogous manner. 

The ratio of the twisted to the untwisted partition function is the vortex free energy 
order parameter. Its behavior as a function of the size of the system characterizes the sys- 
tem's possible phases. By known correlation inequalities it can, furthermore, be related to 
the Wilson and t'Hooft order parameters. Our representations of the twisted and untwisted 
partition functions may now be used to represent the ratio (section [5|). One may exploit 
the parametrization invariance of these representations to ensure that the bulk free energy 
contributions resulting in each decimation step ab^~^ — > ab"^ explicitly cancel between nu- 
merator and denominator in the ratio. One is then left with a representation of the vortex 
free energy solely in terms of an effective action defined on a lattice of spacing a6". 

Now this effective action is constrained by the effective actions corresponding to the upper 
and lower bounds. The latter are easily explicitly computable by straightforward iteration 
of the potential-moving decimation rules. Under successive transformations they flow, for 
space-time dimension d < 4 and any original coupling g deflned at spacing a, to the strong 
coupling regime. This is the regime where the coefficients in the character expansion of the 
exponential of the action become sufficiently small for the strong coupling cluster expansion 
to converge. Confining behavior is the immediate result for the vortex free energy, and, hence, 
'area law' behavior for the Wilson loop (section [6]) . 
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As it is well-known the theory contains only one free parameter, a physical scale which 
is conventionally taken to be (some multiple of) the string tension. This fact comes out in 
a natural way in the context of RG decimations, as we will see in the following. Fixing this 
scale then determines the dependence g{a). The fact that g{a) —> as a ^ is an essentially 
qualitative consequence of the flow exhibited by the decimations. 

Some concluding remarks are given in section [71 



2 Decimations 

We work on a hypercubic lattice A C Z'^ of length in the x'^-direction, = 1, . . . , d, in units 
of the lattice spacing a. Individual bonds, plaquettes, 3-cubes, etc are generically denoted by 
6, p, c, etc. More specific notations such as 6^ or are used to indicate elementary m-cells 
of particular orientation. We use the standard framework and common notations of LGT 
with gauge group G. Group elements are generically denoted by U, and the bond variables 
by Uh € G. In this paper we take G = SU (2). 

We start with some appropriate plaquette action Ap defined on A, which, for definiteness, 
is taken to be the Wilson action 

Ap{Up,P) = I Rexi/2{Up) , Up=l[Ut, (2.1) 

bedp 

with (3 = A/ defining the lattice coupling g. The character expansion of the exponential of 
the plaquette action function is given by 

exp {Ap{U, (3)) = dj Fj{/3) x,{U) (2.2) 
j 

with Fourier coefficients: 

FjiP) = JdU exp{Ap{U,P)) J,Xj{U) ■ (2.3) 

Here dU denotes Haar measure on G, and Xj the character of the j-th representation of 
dimension dj. So, for SU(2), the only case considered explicitly here, all characters are real, 
j = 0, i, 1, |, . . ., and dj = (2j + 1). ([23]) implies that Fq > Fj, ah j / 0. Explicitly, one 
finds 

2 

in terms of the modified Bessel function ly. 

It will be convenient to work in terms of normalized coefficients: 



Fj{(3) = -h^{P) (2.4) 



so that 



exp(Ap([/,/3)) = Fo[l + Y.d,Cj{(5)xm 

^ Fofp{U,f3). (2.6) 
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The (normalized) partition function on lattice A is then 

= fdUAH fpiUp,f3) ^ f dfil, (2.7) 

where dU\ = YibeA dUb, and expectations are computed with the measure df-t\ = dfj,^/Z\{(3). 
The action (j2.ip is such that 

Fj{(3) > , hence 1 > Cj{(3) > all j , (2.8) 

which implies that the measure defined by (j2.7p satisfies reflection positivity (RP) both in 
planes without sites and in planes with sites. Note that lim^^oo Cj (/?) = 1- 

Let A^*^) be the hypercubic lattice of spacing 6"a, with integer b > 2, and Z j^(n) {{cj (n)}) 
denote a partition function of the form (j2.7p defined on A^") in terms of some given set of 
coefficients {cj{n)}: 

Z^(n)({cj(n)}) = j dUf^^n) Y[ [l + ^djCj{n)xj{Up) 

= / ^f^AW n fp{Up,n) = J dfil^^^, (2.9) 
peA(") 

where dUf^(„) = n6eA(") dUb- We also employ the notations 

9p{U, n) = fp{U, n)-l = J2dj Cj{n) Xj{U) , (2.10) 

and II • II for the || • ||oo-iiorm: 

\\g{n)\\=Y^d]cj{n). (2.11) 

One has the simple but basic result: 

II. 1 For Zj^(n){{cj{n)}) given by 112. 9\) with Cj{n) > for all j, and periodic boundary 
conditions, 

(i) Z^{„) ({cj (n)}) is an increasing function of each Cj{n): 

9Z^(„)({Q(n)})/5c,(n) >0; (2.12) 

(u) 

^A(")({c,(n)})> [l + Y,d]c,{nf\ . (2.13) 

(|2.12p is an immediate consequence of RP in planes without sites. The proof of (j2.13p . also 
based on RP, is given in Appendix A. Strict inequality in fact holds in ()2.12p and ()2.13p . with 
equality only in the trivial case where all Cj(n)'s vanish. In particular, one has 

^Aw({c,(n)}) > 1. (2.14) 

Simple as (j2.14p is, it is not trivial, as it requires non-negativity of Cj(n)'s, and will be useful 
in the following. 
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2.1 Construction of decimation transformations 

To perform an RG transformation a — > ba, the lattice is partitioned into li-dimensional 
decimation cells of side length ba. Various approximate decimation transformations may 
be devised involving the 'weakening', i.e. decreasing the Cj's of interior plaquettes, while 
compensating by 'strengthening', i.e. increasing Cj's of boundary plaquettes of each cell. The 
simplest such scheme [5], which is adopted in the following, implements complete removal 
of interior plaquettes. This may be pictured [6] as moving the potentials due to interior 
plaquette interactions to the boundary. 

This 'potential moving' may be performed as the composition of elementary steps. The 
elementary potential moving step is defined in terms of a 3-dimensional cell of side length ba 
in a given decimation direction, say the x'^-direction, and length a in the other two directions 
/i, v. Two such 3-cells adjacent along the K-direction are shown in Figure [TJ The (b — 1) 
interior plaquettes in each cell perpendicular to x'^ (shaded) are removed, i.e. 

Ap{Up) ^ (2.15) 

for the action at their original location, and displaced (arrows) in the positive x** direction to 
the position of the corresponding plaquette (bold) on the cell boundary. There the displaced 
interior plaquettes are combined with the boundary plaquette into one plaquette p with action 
'renormalized' by some appropriate amountH (q: 

Ap{U) ^ Co Ap{U) . (2.16) 

A complete transformation consists of performing this elementary operation successively in 
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Figure 1: Basic plaquette moving operation, 6 = 2 

every lattice direction k = 1, . . . , d in such a way that eventually one is left only with plaquette 
interactions on a lattice of spacing ba. In practice, there is no reason for a choice other than 
6 = 2, but, for clarity, we keep general (integer) 6. The result of a complete transformation 
is given by equations (j2.18p - (j2.22p below, to which a reader may turn directly. 

To describe this process in more detail, let the lattice be partitioned into d-dimensional 
hypercubic decimation cells a'^ of side length ba in each lattice direction. Plaquettes interior 
to a a'^ are defined as those not wholly contained in its {d — l)-dimensional boundary da'^ . 
Consider the effect of successive application of the elementary moving operation to plaquettes 

^One may take this renormalization factor to depend on the move direction, but we need not consider these 
more general transformations here. 
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of fixed orientation, say [fJ-i']- There are {d — 2) normal directions 7^ /i, z^, i = 1, ■ ■ ■ ,d — 2, 
in which a plaquette p^u can be moved. Interior p^,y's in each a'^ are first moved to the cell 
boundary da in groups of (6— 1) parallel plaquettes, along , say, the positive Ki-direction (as 
in Figure 1). They end up in the face cjI'J C da'^ perpendicular to the Ki-axis. There each 
group is identified with the plaquette present at that location and merged in one plaquette 
Pnu G (^ifi with a 'renormalized' action p. 160 ). Similarly, p^^, plaquettes in each face 
o"!'^ C da'^, with i 7^ 1,/U,z/ are moved along the Ki-axis in groups of (6 — 1) to the face 
crlttf ^ C daif^ normal to the ki and Kj directions, where they are merged and renormalized. 

There are now (d — 3) directions inside the {d — l)-dimensional face al'^ in which a 
[/ii/] -plaquette can move. Thus in proceeding to apply the elementary moving operation 
successively in all directions, the once-moved-renormalized p^jy's in ai'^ are next moved, 
in groups of {b — 1) plaquettes in the positive K2-direction, to the face a"i'^K2^ C ^\ 
Similarly, the once-moved-renormalized p^jy's inside a face aifj^nf^ are moved provided K2 is 
among the {d — 4) available directions normal to a [;u;/]-plaquette inside alfj^nf^ . 

Continuing this process in the remaining directions Kj, i = 3, . . . , (d — 2), the set of [^z^]- 
plaquettes on the initial lattice ends up in the 2-dimensional faces (7^^ ^2 K(d_2) '^^kik2 K(rf_3) 

The above process, described for plaquettes of one fixed orientation [f-tu], is carried out 
for each of the d{d — l)/2 possible choices of plaquette orientation [6]. 

The end result of the process is then a lattice having elementary 2-faces of side length 
ba, each tiled by b^ plaquettes of side length a. The action of each of these 6^ plaquettes has 
been renormalized according to ()2.16p by a total factor of 

Cf^ = C . (2.17) 

This is expressed by ()2.22p below. 

The integrations over the bonds interior to each 2-face of side length ba are now carried 
out. This merges the 6^ tiling plaquettes into a single plaquette of side length ba. These 
integrations are exact and do not change the value of the partition function that resulted 
after the completion of the plaquette moving operations. We, however, allow further renor- 
malizing the result of these integrations by introducing, in addition to Co) another parameter, 
r (cf. (j2.19p below). This completes the decimation transformation to a hypercubic lattice 
of spacing spacing ba. 

The important feature of this decimation transformation is that it preserves the original 
one-plaquette form of the action, so the result can again be represented in the form (j2.6p . 
The transformation rule for successive decimations 



a ^ba ^ b^a^ > b'^'^a 6" 



a 



A ^ a(^) ^ A(2) ^ . . . ^ A(-i) ^ A(") ^ • • • , 

is then: 

fp{U,n- 1) ^ Fo(n) UU,n) = Fo{n) [l + ^ dj c,{n) Xj{U)] . (2.18) 
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Cj[n) 



Cj{n) 
Fo(n) = Fo{n) 



fe2 



The n-th step coefficients FQ{n), Cj{n) are obtained from the coefficients Cj{n — 1) of the 
previous step by 

(2.19) 
(2.20) 

(2.21) 
(2.22) 



where 
and 



c,in) = Fjin)/Fo{n)<l, j^O, 



dU 



/(C/,n-l)]^lxi(^). 



The n = coefficients are the coefficients Cj{(3) on the original lattice A. ()2.22p encodes the 
end result of the plaquette moving - renormalization operations described above, with C, of 
the form (j2.17p : and (j2.19p . (j2.20p that of the subsequent 2-dimensional integrations, and 
further renormahzation by the parameter r. 

It is easily seen that fp{U,n) > given that this holds forn = (cf. <^M)- The 

effective plaquette action on lattice A^") of spacing b'^a is then 

fp{U,n) = [l + J2djc,{n)x,{U)] (2.23) 

= exp(^p(C/,n)) , (2.24) 
with effective couplings defined by the character expansion 

Ap{U,n) = (3o{7i) + J2Pi{n)d^Xi{U) . (2.25) 

A point on notation. In the above we used the notations FQ{n), Cj{n), f3j{n), etc, which 
do not display the full set of explicit or implicit dependences of these quantities. Thus, a 
more complete notation is: 

Cj{n) = Cj{n,b,C,r,{cj{n- I)}) 

Fo{n) = Fo{n,b,C,{cj{n-l)}) (2.26) 

Dependence on the original coupling /? comes, of course, iteratively through the coefficients 
{cj (n — 1)} of the preceding step. Because of the iterative nature of many of the arguments in 
this paper several explicit and implicit dependences propagate to most of the quantities used 
in the following. To prevent notation from getting out of hand we generally employ short- 
hand notations such as those on the l.h.s. of (I2.26p . unless specific reference to particular 
dependences is required. 

The resulting partition function after n such decimation steps is: 



ZA{P,n) = n Fo(m)l^l/^'"' ZA(„)({c,(n)}) , 



(2.27) 



771=1 



with Z ^(n) {{cj (n)}) of the form (|2.9p and coefficients (j2.26p resulting after n steps according 
to (|2.19p - (|2.22p . The bulk free energy density resulting from decimating from scale a to 
6"a is then ^J^^^ lnFo("T-)/^'"°' ) each term in this sum representing the contribution from 
^(771-1) _^ ^777^ spgcifigfi by (j2.18p . The partition function (j2.27p is, of course, not equal 
to the original partition function Z\{(3) of (12. 7p since the decimation transformation is not 
exact. How they are related will be addressed below. 
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2.2 Some properties of the decimation transformations 



The transformation rule specified by (j2.18p - ()2.22p is meaningful for real positive Here, 
however, a basic distinction can be made. As it is clear from (j2.22p . for integer the important 
property of positivity of the Fourier coefficients in (j2.18p is maintained at each decimation 
step: 

Fo{n)>l, l>Cj(n)>0 (integer C) • (2.28) 

This means that reflection positivity is maintained at each decimation step. This clearly is 
not guaranteed to be the case for non-integer C. Thus non-integer ^ results in transformations 
that, in general, violate the reflection positivity of the theory (assuming a reflection positive 
original action). 

It is important in this connection that, after each decimation step, the resulting action 
retains the original one-plaquette form, but will generally contain aZ/ representations in (j2.25p . 
Furthermore, among the effective couplings (3j{m) negative ones will occur. These features 
are present in general, even after a single decimation step a — > 6a starting, as we did, with 
the single (fundamental) representation Wilson action ()2.ip . For integer however, the 
resulting effective action ()2.25p . even in the presence of some negative couplings, still deflnes 
a reflection positive measure, since, as just noted, the expansion of its exponential (j2.23p 
gives positive coefficients (12.281) . 

It is also worth noting that, given a set of initial coefficients (j2.5p . the transformation rule 
(j2.19p - (j2.22p with integer Q can be explicitly evaluated, to any desired accuracy, by purely 
algebraic operations, namely repeated application of the KG reduction rule 

i+j 

XiXj= E Xk (2.29) 

k=\i-i\ 

in ()2.22p and character orthogonality - no actual integrations need be carried out. 
The choice (cf. dHZD) 

Co = 1 + (6-1) =^ C = b'-' (2.30) 

is special. It increases the couplings of receiving plaquettes, at each basic moving step, by an 
amount exactly equal to that of the corresponding displaced plaquettes. This, together with 
r = 1, is essentially the original choice in [5] as reformulated in [6], and will be referred to as 
MK decimation. It will be important in the following!! 

There are various other interesting properties of the decimations that can be derived 
from (j2.19p - (j2.22p . The following one is particularly important. The norm (j2.1ip of the 
coefficients obtained by application of (I2.19|) - (I2.22P with integer C ^ 1 ^'^d r = 1 satisfies 
(Appendix D): 

{C-l)62 



\\g{n + l) 
Assume now that 







b2 




< 


\CMn)\\ 




"l + lb(n)ir 



2.31 



|5(n)|| <exp(-C„). a>0, (2.32) 



^It is worth noting in this context that in numerical investigations of the standard MK recursions in gauge 
theories [8] fractional 6, (1 < 6 < 2), which by (|2.30p corresponds to non-integer has often been used. 
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fe2 



for some n. Then 

\\g{n + 1)11 < C\\g{n)\\ exp 
where = In C + (C — 1). The recursion 

Cn+l 

gives 

Cn ~ 



((C-i)6^ 



< 



exp 



iCn-k)b^ 



Cnb" 
b^k 



62 - 1 



b^k 



62 - 1 



Cn+m. 

II. 2 If for some n the norm of coefficients \2.11\) obeys i2. 3^) with 

b^k 



Cn> 



62 - 1 



\g{n + m) 



(2.33) 

(2.34) 
(2.35) 

(2.36) 

-> as 



then, under iteration of the decimation transformation 12.19\) - 
m ^ oo according to ^2.35\) . 

This fah-off behavior is immediately recognizable as "area-law" . 

If one assumes that Cj{n) are small enough so that the theory is within the strong coupling 
regime, this behavior can be immediately deduced for the leading coefficient Ci/2(f^) directly 
from (I2l9l) - ([2:22]) : 



ci/2(n + 1) = c^/2{nf exp ( [InC + 0(ci/2(n)) ] b^ 



(2.37) 



The result (I2.33P gives then an estimate of the corrections due to all higher representations. 
What is noteworthy here, however, is that the condition (|2.36p is rather weaker than the 
commonly stated conditions for being inside the convergence radius of the strong coupling 
cluster expansion (cf. section [6]). 

We note two further properties of the decimation transformations (I2.18P - ()2.22p . The 
first is that with r = 1 they become exact in space-time dimension d = 2 since then, from 
(j2.17p . C = 1- The second is that, with C, = b^'^~'^\ vanishing coupling (7 = is a fixed point in 
any d, i.e. MK decimation is exact at zero coupling. This follows simply from the fact that 



lim 



ess. sup e 



of\ 



for any normalized measure dv[x). Applying this to the result of performing the plaquette 
moving operation starting from (|2.6p . and with p' £ A labeling the plaquettes tiling the 
plaquettes p G A^^^ , one has 



lim 

/3-*oo 



n 

pGA(i) 



jdUA n ne^p(/5^^'"'^^xi/2(f/p 



exp 



lim 

/3-»oo 



dU\ Y[ exp (/3 ^Xi/2(t^p 



peA 



This clearly holds also for r 7^ 1, as is evident from the fact that lim^_^oo Cj(/3) 
fixed point is easily seen to be unstable. 



.(2.38) 
1. This 
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3 Partition function 



Since our decimations are not exact RG transformations, the partition function does not in 
general remain invariant under them. The subsequent development hinges on the following 
two basic propositions that relate partition functions under such a decimation. 

3.1 Upper and lower bounds 

Consider a partition Z^(„-i) on lattice A^""^) of the form (j2.9|) given in terms of some set of 
coefficients {cj{n — 1)}. Apply a decimation transformation (j2.18p - (j2.22p performed with 
= b^'^~'^\ Denote the resulting coefficients by c'j , , i.e. 

cf(n,r) EE c,(n,6, C = fo('-'\r,{c,(n-l)}) (3.1) 

F^{n) ^ Fo(n,6, C = 6('^-'),{c,(n-l)}). (3.2) 

Note that 

c'^{n,r) = cf{n,lY. (3.3) 

111.1 For Z^(n~i) of the form /i2.9\) . a decimation transformation 112.18]) - 112. 22\) with 
= b'^~^ and < r < 1 results in an upper bound on Zj^(„-i) : 

^A("-)({c,(n-l)}) < F^{n)\^'"'\Z^,.,{{cf{n,r)}). (3.4) 

The r.h.s. in \3.4\ ) is a monotonically decreasing function of r on < r < 1. 

Given partition function Zy^(„-i) on lattice A^"'"^) of the form (|2.9p in terms of some set 
of coefficients {cj{n — 1)}, let 

c^(n) = Cj{n-lf (3.5) 
F^{n) ^ 1. (3.6) 

111.2 For Zj^((„^iy) , Zj^(n) of the form \2. 9\) : 

^Aw({c-'(n)}) < Z^(„-,({c,(7^ - 1)}) . (3.7) 

The proof of III.l is given in Appendix A, where somewhat stronger results than ()3.4p 
are actually obtained. III. 2 is a corollary of (j2.13p (Appendix A). For the argument in the 
rest of this paper, the precise form of the lower bound is in fact not important. By Il.l(i) a 
further lower bound is obtained by replacing cj'(n) in III. 2 by, for example, 

c^(n) = Cj{n - lfc^j{n,r) (3.8) 

since < Sj{n,r) < 1. Another choice is to simply set 

c|(n) = 0, (3.9) 
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which is a restatement of (|2.14p . 

A related lower bound, which, in analogy to the upper bound in III.l, can be formulated 
directly in terms of the transformations ()2.18p - p.22p . is obtained by taking Cj'{n) in III. 2 
to be given by: 

Cj{n) = Cj{n,b, C =l,r =l,{cj{n-l)}) 

= Cj{n-1)^\ (3.10) 
Fo^(n) = Fo(n,6, C = l,{cj(n-1)}) 

= 1. (3.11) 

With this choice of Cj{n), note that III.l - III. 2 imply the fact that the decimations (|2.18p - 
(j2.2ip become exact for d = 2 and r = 1. III.l says that, after removal of interior plaquettes, 
modifying the couplings of the remaining plaquettes by taking C, = 6"^"^ (and r < 1) results 
into overcompensation. III. 2 says that decimating plaquettes while leaving the couplings of 
the remaining plaquettes unaffected (C = i" = 1) results in undercompensation. The proof 
of III.2 for cj^(n) given by (I3.10|) is similar to that of II. 1, but need not be given here, since 
the weaker bounds above will suffice. 

In the following it will in fact be more convenient to take ()3.8p or (j3.9p for the definition of 
the lower bound coefficients c^(m). Use of the stronger lower bounds above may be preferable 
for numerical investigations, but does not contribute anything further to the argument in this 
paper. 

III.l and III. 2 give upper and lower bounds on the partition function after a decimation 
step. It is then natural to interpolate between these bounds. 



3.2 Interpolation between upper and lower bounds 

Introducing a parameter a G [0, 1], we define coefficients Cj{m, a, r) interpolating between cj" 
at Q = and (j3.ip at a = 1: 

Cj{m,a,r) = {1 — w{a)) Cj (m) + w{a) {m,r) , < r < 1. (3-12) 

with 

w{0) = , w{l) = 1 , w'{a) > . (3.13) 

For example, 

w{a) = (3.14) 
e — 1 

There is clearly a variety of other choices than (|3.12p for these interpolating coefficients. We 
always require that 

dcj{m,a,r)/d a> Q , (3.15) 

which is satisfied by (I3l2]) - (I3l3]l . 

Similarly, we define coefficients interpolating between (13. 6p and p.2p . For our purposes 
it will be convenient to take 

Fo(m, h, a, t) = F^imf'^"^ , (3.16) 
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where ht denote a family of monotonically increasing smooth functions of a, labeled by a 
parameter t £ [ta, tb] , and such that 

ht{0) = 0, ht{l) = l. (3.17) 

We write ht{a) = h(a,t). Examples ar^ 

h(a,t) = ex.p(—a(t) ^ , h(a,t) = a'^^^\ h(a, t) = tanh ( -) , 

\ a J a[t) (1 — a) 

0<a<l, <ta <t <tb < oo, (3.18) 

where a{t) is a smooth monotonically increasing positive function on [iai^ft]) e.g. cr{t) = t. 
The interpolating partition function on A^*") constructed from cj and Fq is now defined 

by 

Zj^im){P,h,a,t,r) = Fo(m, /i,a,i)l'^''"'l Z^(™)({cj(m,a,r)}) (3.19) 

where 

= /'^^A(-) n fp{Up.rn,a,r). (3.20) 

Combining II. 1, (|3.15p and the fact that Fq is, by definition, also an increasing function of a 
one has 

III. 3 The interpolating free energies In Z^(m){{cj{m, a, r)^) and \nZ^{m){l3,h,a,t,r) are 
increasing functions of a: 

ainZ^t™, ({cj(m, a, r)}) /5a > 0. (3.21) 



Equality in ()3.2ip applies only in the trivial case were all the coefficients c^-'s vanish. 
In terms of (l3lB . III.l and III.2 give 

(/3, h, 0, t, r) < ^^(m-i) < (/3, h, 1, t, r) . (3.22) 

Now Zf^i^m) (/3, /i, a, t, r) is continuous in a. It follows from (j3.22p that there exist a value of 
a in (0, 1): 

a(m, /i, t, r, {cj (m — 1)}, 6, A) = a^j] (*> (3.23) 

such that 

^^(m)(^,/i,ai"2(*''^)'*''^) = -^A('"-i) • (3-24) 

In other words, at each given value of t, r, there exist a value of a at which the partition 
function on A^"*), resulting from a decimation transformation A^'""^) — > K^^\ equals the 

^Supplementing these definitions at a = as needed is understood. Thus, h(oi, t) = on q < in the first 
example in p.l8|l : and standard smoothing in the second example: replace a in /i by g,,{Q) — J p^{a — x)g{x)dx, 
where g{x) — x for a; > 0, g{x) — for x < 0, and Pt{x) is C°°, has support inside \x\'^ < and satisfies 
pe > and J pe — 1. 
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partition function on A^"^ ^\ This value is unique by III. 3. By construction, a^(f,r) is 
such that ()3.24p remains invariant under variation of t, r in their domain of definition, i.e. 
a^^(i,r) represents the level surface of the function Z^(m){f5,h,a,t,r) fixed by the value 
Z^(m-i) ■ The parametrization invariance under varying t will be important later. 

We now examine the dependence on t, r in (I3.23|) more closely. Given Z^{m-i) and some 
interpolation h, assume that (I3.24|) is satisfied at the point (to,ro,a = Oaj])- Then, by the 
implicit function theorem, applicable by III. 3, there is a function a^^(t,r) with continuous 
derivatives such that a"^'j^{tQ,rQ) = aj^^' uniquely satisfies (13.24p in a sufficiently small 

neighborhood of (to, ^'Oj aj^^)- But since a solution to (|3.24p exists for each choice of t,r 
in their domain of definition, this neighborhood can be extended by a standard continuity 
argument to all points of this domain. a^^(t,r) then represents the regular level surface of 
the function (|3.19p fixed by (|3.24p . Furthermore, 



dt 



v{atl{t,r),t,r) , (3.25) 



where 



with 



— h ^A(m) (a, r) 

^A(™,(a,r) = ^-^^p^P^— lnZ^(,„) ({c,(m, a, r)}) >0. (3.27) 

We will always assume that h is chosen such that dh/dt is negative. This is the case with 
the examples (|3.18p . Then, from (j3.26p . f>OonO<a<l, with u = at a = and q = 1. 

It is also useful to equivalently view ot^xh{t, r) as the solution to the ODE 



da/dt = v{a,t,r), a e (0, 1) , (3.28) 

(m) 



a{to) = al")] > , tQ£[ta,tb]. 



Then standard results of ODE theory imply the existence of a unique solution in a neighbor- 
hood of a]^ h > 0' which can in fact be extended indefinitely forward for all t > too 
A short computation using ()3.25p gives 

d/i(a^(t,r),t) da'i';lit,r) 

-Jt = ^A('^)io'A,hit,r),r), (3.29) 

as it should for consistency with p.24p . (j3.29p and (j3.25p make apparent what the effect of a 
parametrization change due to a shift in t is. Increasing (decreasing) t increases (decreases) 
the contribution of In Zj^(m){{cj {in, a^^l{t,r),r)}) while decreasing (increasing) by an equal 

amount the contribution from lnFg^(?Ti) /i(q^ (t, r), t) lA^^^I, so that the sum stays constant 
and equal to InZ^(m-i) in accordance with ()3.24p . 



^Indeed, v is differentiable on a^^' < a < 1 and vanishes at a 
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The derivative w.r.t. r is similarly given by (jB.ip - ()B.2p in Appendix B. Now, by (III.l), 
the upper bound in (I3.22P is optimized for r = 1, which would appear to make consideration 
of other r values unnecessary. The reason one may want, however, to vary r away from unity 
is the following. 

The values a^^l^t, r) lie in the interval (0, 1). Consider the possibility that one finds that 

a^^(tm;l) differs from 1 only by terms that vanish as the lattice size grows. This means 

that, since > in (j3.25p . aj("^^(t, 1) is, to within such terms, a constant function of t for all 
t > tm.. For the purposes of the argument in the following sections we want to exclude this 
possibility, and ensure that, at least in some neighborhood of a chosen t value, the derivative 
(j3.25p is non-vanishing by an amount independent of lattice size. 

We require that 

y < aH(t,r) < 1-5 , (3.30) 

with 6 > 0, 6' > independent of the lattice size |A(™)|. The lower bound requirement is 
easily shown (Appendix B) to be automatically satisfied by combining II. 1 and (I3.24p . As it 
is also shown in Appendix B, one may always ensure that the upper bound requirement in 
(j3.30p holds by choosing the decimation parameter r to vary, if necessary, away from unity 
in the domain 

l>r>l-e, (3.31) 
where < e <C 1 with e independent of |A(™)|. 

With (13301) in place, (13:25]) and (13:29]) imply (Appendix B) that 

^(t,r) > vi{S) > 0, -^(ay(t,r),t) > m{S) > 0, (3.32) 

where r]i, r]2 are lattice-size independent. Furthermore, if (|3.30p already holds for r = 1, it 
also holds for any r in (|3.3ip . We may as well then simplify matters in the following by setting 
the parameter r to the value r = 1 — e with some fixed small e. This e may eventually be taken 
as small as one pleases after a sufficiently large number of decimations have been performed. 
This has an obvious meaning in the context of iterating the decimation transformation as 
pointed out in subsection 13.41 below. We accordingly simplify notation by dropping explicit 
reference to r, except on occasions when a statement is made for general r values. Thus we 
write Cij^\{t) = Cij^\{t, 1 — e), c^{m) = c^{m, 1 — e)), etc. 



3.3 Representation of the partition function on decimated lattices 

So, starting on the original lattice spacing a, with partition function given in terms of co- 
efficients {cj{f3)}, one may iterate the procedure represented by (|3.22p - (j3.24p . Taking the 
same interpolation family h in every cycle, an iteration cycle consists of the following steps. 

(i) A decimation transformation A^™"^) A^™) given by the rules (|2.18p - (|2.22p applied 
to the coefficients in Zj^(m-i), and resulting into the upper bound coefficients on A^™-* 
according to (|3.ip - (|3.2p and (|3.4p . Similarly, a lower bound on A^*") is obtained 
according to ()3.7p with lower bound coefficients given by ()3.6p and ()3.8p or ()3.9p . 
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(ii) Interpolation between the resulting upper and lower bound partition functions on A^*") 
according to (f3l2]) . (fCTIl . and (f3l9]) . (fCTj) . 

(iii) Fixing the value < oi'Ahit) < 1, eq. (|3.23p . so that the (m — l)-th step partition 
function Zj^^(m~i) is preserved, eq. ()3.24p . 

(iv) Picking a value of the parameter t = tm, to fix the coefficients {cj{m, aj^^(tm))} of the 
resulting partition function Zj^(,n), and return to step (i). 

This scheme for the coefficients in ^^(m) may be depicted as follows: 

cm 



< 



< 



{c'/{2)} 



(3.33) 



The result after n iterations is then: 
Za(/3) = fdUAll fp{U,(3) 



|A|/b" 



m=l 



(3.34) 



ZA(„)({c,(n,ai';Uin))})- (3-35) 



(j3.35p is an exact integral representation on the decimated lattice A^") of the partition function 
Za originally defined on the undecimated lattice A by the integral representation (12. 7p or 

III. 3 allows the iterative procedure leading to ()3.35p to be implemented in a slightly 
different manner, one that turns out later to be more convenient for our purposes. Since by 
III.3 



^A(-) ({cj{m,a^J^l{tm))}) < ^AC-) ({cj("i, 1)}) = Zac^) ({cf (m)} 



(3.36) 



an upper bound for each successive iteration step is also obtained by applying III.l to the 
r.h.s. rather than the l.h.s. of the inequality sign in (|3.36p . The only resulting modification 
in the above procedure is in step (i): the upper bound coefficients c^(m) and Fq {m) on A^*") 
are computed according to (j3.ip and (|3.2p but now using the set {c^ {m — 1)} rather than the 

set {cj{m — 1, Oa™^ ^\tm-i))} as the coefficient set of the previous step. 

The same alternative can be applied to the lower bounds in (|3.33p . Since, again by III. 3, 
one has 



^Aim) [{cf{m)}^ = Z^^m) ({Cj(77T,,0)}) < [{Cj {771, J^litrr,))}) , 



(3.37) 
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a lower bound for each successive iteration step is also obtained by applying III. 2 to the l.h.s. 
rather than the r.h.s. of the inequality sign in (j3.37p . If one adopts p.9p . this makes no 
difference since the lower bound coefficients equal zero at every step. If one uses (j3.8p . the 
resulting modification to (j3.33p is that in step (i) the lower bound coefficients c^(m) on A^*"-* 

are now computed using the set {cj'(m — 1)} rather than {cj{m — l,a^^f^ ^\tm~i))} as the 
coefficient set of the previous step. 

One may adopt either or both modifications following from p.36p or (I3.37p . Adopting 
both, the iterative scheme for the coefficients in Z^{m) replacing ()3.33p is: 




1 I (3.38) 

{c[(2)} < {c,(2,af < {cf(2)} 



This again leads, after n iterations, to the representation ()3.35p . Note, however, that the 
actual numerical value of a^^l^{tm) in (j3.35p . fixed at each step by requiring (j3.24p . will, in 
general, be different depending on whether scheme (j3.33p or (j3.38p is used for the iteration. 
Also note that the upper bounds (m) in (j3.38p are not optimal compared to those in 
(|3.33p . The scheme (|3.38p . however, turns out to be more convenient for our purposes in the 
following. 

3.4 Discussion of the representation ( 13.351) 

As indicated by the notation, on any finite lattice, the values possess a lattice size 

dependence. This weak dependence enters as a correction that vanishes inversely with lattice 
size. Indeed, by the standard results on the existence of the thermodynamic limit of lattice 
systems, for a partition function Z^{m)({cj}) of the form ()2.9p on lattice A^™) with torus 
topology (periodic boundary conditions): 

lnZA(,„)({c,}) = lA^'")] (^({c,}) + 6^j,^^,{{cj}) , (3.39) 

f{{cj}) being the free energy per unit volume in the infinite volume limit, and 5ipj^{m) {{cj}) < 
0(constant)|£| From this and (|3.24p it is straightforward to show that 

it, r) = ) (i , r) + Ja^ (t , r ) (3.40) 

®That is, there are no 'surface terms' for torus topology. In fact surface terms arising with other, e.g. free, 
boundary conditions can be precisely defined as the difference in the free energies computed with periodic 
versus such other boundary conditions [11] , 
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with 6a'x l{t,r) as some inverse power of lattice size in the large volume limit. In fact, we 

have already established the presence of a lattice-size independent contribution in aj[^(t,r) 
in an alternative manner through (j3.30p . i.e. the fact that in (j3.40p one must have 

4'"^(t,r)>y. (3.41) 

An explicit expression for 6' is given by (|B.17p . (|B.19p . 

At weak and strong coupling the {t, r) values may be estimated analytically by com- 
parison with the weak and strong coupling expansions, respectively. In general, starting from 
(I3:22D . the location of J'^l satisfying ^^TM may be formulated as the fixed point of a con- 
traction mapping. This allows in principle its numerical determination, for given values of 
all other parameters, to any desired accuracy. 

For our purposes here, however, the actual numerical values of the a^^'s, beyond the fact 
that they are fixed between and 1, will not be directly relevant. The main application of the 
representation (j3.35p in this paper will be to relate the behavior of the exact theory to that 
of the easily computable approximate decimations bounding it without explicit knowledge of 
the actual values. 

It is important to be clear about the meaning of (j3.35p . The partition function Z\{f3) 
is originally given by its integral representation (13.34p on lattice A of spacing a. (13.35P 
then gives another integral representation of Z\((3) in terms of an integrand defined on the 
coarser lattice A(") of spacing b"'a plus a total bulk free energy contribution resulting from 
decimating between scales a and 6"a. The action Ap{U,n,a^^\) in ({c^ (n, a^l^}) is 
constructed to reproduce this one physical quantity, i.e. the free energy In Z/^{(3), nothing 
more and nothing less. In particular, it is not implied that this action on A^") can also be used 
to exactly compute any other observable. For that one would need to attempt the previous 
development from scratch with the corresponding operator inserted in the integrand. 

Recall that, by (I3.15p . the coefficients Cj(n, a, r)'s are increasing in a, and Cj{n,l,r) = 
c^{n,r), Cj{n,0,r) = Cj{n): 

cf{n) < c,{n,ai:l{t))< {n) , < ai'^f (t) < 1 . (3.42) 

Thus, the coefficients Cj{n,a'^\{t) ) in the representation ()3.35p are bounded from above by 

c^(n) no matter what the actual values of a^|j(i) are. 

When considering the implications of this bound under successive decimations the ad- 
vantage of employing scheme (|3.38p . rather than (j3.33p . becomes clear. The coefficients 
(n) on the r.h.s. column in (13.38P are obtained by straightforward iteration of the deci- 
mation rules ()2.19p - (j2.22p with ( = 6*^"^; i.e. only knowledge of the c'j{n — 1), not of the 

Cj(n— 1, a^^f^ ^\tn~i)), is required to obtain the c^(n) at the n-th step. The flow of these {n) 

coefficients then constrains the flow of the exact representation coefficients Cj(n, oi^i^\{tn)) ac- 
cording to (j3.42p from above. In particular, if the (^{n)'s on the r.h.s. column in VJ.'J^) 
approach the strong coupling fixed point, i.e. 

F^{n) 1, (n) ^ 0, as n ^ oo , (3.43) 
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so must the Cj{n,a^'fJ 's in the representation 113. 35\) V\ 

Now the coefficients c'j {n, r) at r = 1 are the MK decimation coefficients (cf. section [2^ . 
As it is weh-known, the MK decimations for SU{2) (and also SU{3)) are found by exphcit 
evaluation to indeed flow to the strong coupling fixed point ()3.43p for all starting /3 < oo and 
d < 4. Above the critical dimension d = 4, the decimations result in free spin wave behavior 
{c^ {n, 1) — > 1 as n ^ oo) starting from any /3 > /3q, where f3o = 0(1). 

Here, for reasons discussed at the end of section [3121 we take r in the range (|3.31|) . This 
may be viewed as fixing the direction from which the point ( = b^^~'^^ , r = 1 in the parameter 
space of the iteration (12.19P - (|2.22p is approached. This is actually irrelevant for the flow 
behavior of the {n, 1 — e) = cj {n) since, in the case of SU (2) considered here, this point is 
a structurally stable point of the iteration^ 

Note that zero lattice coupling, g = 0, is a fixed point as it is for the MK decimations. 
This is also evident from lim^^oo Cj(/3) = 1 and III. 2. 

What does (j3.43p combined with (j3.42p imply about the question of confinement in the 
exact theory? The fact that the long distance part, Z^(n)({cj(n, a^|^)}), in (j3.35p flows in 
the strong coupling regime does not suffice to answer the question. It is the combined contri- 
butions from all scales between a and b'^a in (j3.35p that add up to give the exact free energy 
IuZa^P). Indeed, recall that, by a parametrization change by shifts in t at each decimation 
step, one can shift the relative amounts assigned to these various contributions keeping the 
total sum flxed (cf. remarks immediately following (j3.29p . This parametrization freedom will 
in fact be important in the following. On the other hand, the fact that by (j3.42p the flow of 
Cj{n, 0(^^\{tn)) to the strong coupling regime is independent of such parametrization changes 
is strongly suggestive. At any rate, to unambiguously determine the long distance behavior 
of the theory one needs to consider appropriate long distance order parameters. 



4 'Twisted' partition function 

The above derivation leading to the representation (j3.35p for the partition function cannot be 
applied in the presence of observables without modification. Thus, in the presence of operators 
involving external sources, such as the Wilson or 't Hooft loop, translation invariance is lost. 
Reflection positivity is also reduced to hold only in the plane bisecting a rectangular loop. 
Fortunately, there are other order parameters that can characterize the possible phases of 
the theory while avoiding most of these complications. They are the well-known vortex free 
energy, and its transform with respect to the center of the gauge group (electric flux free 
energy). They are in fact the natural order parameters in the present context since they 
are constructed out of partition functions, i.e. partition functions in the presence of external 
fluxes. 

Let Z\{t^u,(3) denote the partition function with action modifled by the 'twist' t^^, i-e. 
an element of the group center, for every plaquette on a coclosed set of plaquettes V^u winding 

^ To strictly draw the same conclusion from the alternative scheme (|3.33|) requires an additional step, such 
as showing that the Cj(n)'s computed according to the scheme (|3.33|) flow to the strong coupling regime if 
those computed according to (|3.38|) do. 

*It is, however, very much relevant in cases where this point is not structurally stable, e.g. in (7(1). 
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through the periodic lattice in the {d — 2) directions perpendicular to the fi, and z^-directions, 
i.e. winding through every [/xi/]-plane for fixed fj,, u: 



if p£V, 



(4.1) 



A nontrivial twist (t^^, ^ 1) represents a discontinuous gauge transformation on the set V^jy 
with multivaluedness in the group center. Thus, for group SU{N), it introduces vortex flux 
characterized by elements of 'Ki{SU {N) / Z {N)) = Z{N). The vortex is rendered topologically 
stable by being wrapped around the lattice torus. 

In the case of SU(2) explicitly considered here, there is only one nontrivial element, 
T^u = —1- As indicated by the notation Z\{t^ui P), the twisted partition function depends 
only on the directions in which Vfj,u winds through the lattice, not the exact shape or location 
of V^j/. This expresses the mod 2 conservation of flux. Indeed, a twist t^^, = — 1 on the 



plaquettes forming a coclosed set V^j/ can be moved to the plaquettes forming any other 
homologous coclosed set V ' by the change of variables Ub —Uh for each bond 6 in a set 
of bonds cobounded by V U V ', leaving Z\{t^u,P) invariant. By the same token, Z\{T^u,f3) 
is invariant under changes mod 2 in the number of homologous coclosed sets in A carrying a 
twist. In the following, for deflniteness, we flx, say, ^ = 1, u = 2, and drop further explicit 
reference to the fi, v indices. Also, we write Z\{—l,P) = 

(j4.ip implies that ^ is obtained from Za by the replacement 



fp{Up,a) fp{-Up,a) 



l + Y^[-lf^d,c0)xjiUp 



for each peV , (4.2) 



in (j2.7p . (12. 6p . i.e. only half-integer representations on plaquettes in V are affected. In general 
then, the twisted version of the partition function (j2.9p on A^") is 



(4.3) 



with 



(4.4) 



f^-\Up,n) = [l + ^(-l)2.-^.[V]rf^.,^.(^)^^.(f/^ 

Sp\y] denotes the characteristic function of the plaquette set V, i.e. Sp\y] = 1 if p G V, and 
^p\y] ~ otherwise. A simple result (Appendix A) of obvious physical signiflcance is: 

IV. 1 With Cj{n) > 0, all j, 

Z^~i,{{cj{n)}) < Z^Mc,{n)}) . (4.5) 

Strict inequality holds in fact in (j4.5p for any nonvanishing j3 on any flnite lattice. 

Application of the decimation operation deflned in section[2]on some given Z^(2_i) of the 



form (14.31) results in the rule 



/(-)([/,m-l) 



i^o(m) /(-)([/, m) 



Fo(m) 



1 + 



-l)2i^^M(i,c,(m)x,(C/)l, (4.6) 
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with coefficients Fo(m), Cj{m) computed according to tlie rules (|2.19|) - (j2.22p . Starting on 
lattice A, the twisted partition function resulting after n such steps is 

n 

4"^(/?'") = n Fo{m)\^\/'"'' Z^^,l{{c,{n)}) . (4.7) 

m=l 

Note that the flux is carried entirely in -^^(i)- Indeed, bulk free energy contributions from 
each A(™~^) — > A^™) decimation step arise from local moving-integration operations within 
cells of side length b on A^'^~^\ i.e. topologically trivial subsets, and are thus insensitive to 
the flux presence. The evolution with n of the effective action in Z^^l^ then determines the 
manner in which flux spreads, which is characteristic of the phase the system is in. 



4.1 Upper and lower bounds 

In the presence of the flux, the measure in (j4.3p possesses the property of reflection positivity 
only in (d— l)-dimensional planes perpendicular to any one of the directions p 7^ 1, 2 in which 
V winds around the lattice. One way of dealing with this is to simply consider the quantity 

^A(n)({c,(n)}) ^ i(ZA(„)({c,(n)}) +4;),({c,(n)})) (4.8) 

instead of -^^(i) • It is indeed easily checked that reflection positivity holds for the measure 
in -^^(,1) in all planes. A direct consequence of this (Appendix A) is then the analog of II. 1: 

IV. 2 For •^^(„)({cj(n)}) given by |^.<g| j with Cj{n) > for all j, and periodic boundary 
conditions, 

(i) -^^(„) ({cj(n)}) is an increasing function of each Cj{m): 

5Z+„)({Q(n)})/ac,(n)>0; (4.9) 

(n) 

Ztai^jin)}) > [l + Y.d]c,inf] . (4.10) 

Again, in these bounds equality holds only in the trivial case where all Cj(n)'s vanish. In 
particular, one has 

^i„)({c,(n)}) > 1. (4.11) 

Note that these bounds are identical to those in II. 1. This signifies the obvious fact that 
they bound from below by underestimating the bulk free energies proportional to the lat- 
tice volume, whereas the lattice size dependence of the free energy discrepancy between 
Z^(„) ({cj(n)}) and ^(i) ({cj (n)}) is much weaker. 

Upper and lower bound statements analogous to III.l and III. 2 can be obtained for Zjl^^^y 
One has: 

IV. 3 For of the form ^.8\ ), a decimation transformation ^.6^ , i2.19\) - \2.2S\) with 

( = b'^^^ and < r < 1 results in an upper bound on .• 

^+„_.)({c,(n-l)}) < Fo^(n)l^'"'lz+„)({c7(n,r)}) . (4.12) 
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The r.h.s. in ^.IB^ is a monotonically decreasing function of r on < r < 1. 
IV.4 For of the form g^j; 

Zl,„,{{c'-{n)}) < Z+„_,({c,(n - 1)}) , (4.13) 

where the coefficients Cj{n) are given by 113. 5\) . 

The proof of IV. 3, as well as that of IV.4, an easy corollary of IV. 2, are given in Appendix 
A. It then follows from (|4.9|) ) that (|4.13|) holds also with coefficients Cj{n) given by (|3.8|) or 
(j3.9|) . Again, in analogy to III. 2, IV.4 also holds with cj" given by (|3.1(J|) . but this form will 
not be used here. 

4.2 Representation of Z\ + on decimated lattices 

The procedure of section [3] leading to the representation (j3.35p for Za can now be applied to 
= {Z\ + Zjy ^)/2. One introduces the interpolating coefficients Cj{m,a,r) given by eq. 
(|3.12|) . and Fq^tu, h, a, t) given by eq. (|3.16|) for some choice of interpolation function h such 
as given by the examples (j3.18p . The quantity corresponding to (j3.19p is then given by 

/i, a, t, r) =FoK^, I ^+(™)({c,(m, a, r)}) (4.14) 

where 

Zt(n.){{cj{m,a,r)}) = ^ (z^f™, ({cj(m, a, r)}) + Z^ji ({^j (m, a, r)})) (4.15) 

with Z^(,„)({cj(m,a,r)}) given by ([3:20]) and ({cj(m, a, r)}) given by (jM]) - with 

coefficients Cj{m,a,r). We then have the analog of III. 3: 

IV. 5 The interpolating free energies In Z'^^^,^^^{{cj{m, a, r)}) and In Z^^^j (/3, /i, a, t, r) are 
increasing functions of a: 

5 In Z+„)({c,(m, a, r)}) /5a > 0. (4.16) 

In terms of IV. 3 and IV.4 give 

Z+„)(/3,/i,0,t,r) < Z+_,, < Z+^,(/3,/i,l,t,r). (4.17) 
which implies that there exist a value of a in (0, 1): 

a'^{m,h,t,r,{cjim - l)},fe,A) = a^'^hht^r) (4.18) 

such that 

Z+„)(/3,/i,a+y (t,r),t,r) = Z+_„ . (4.19) 

This value is unique, for given values of t, r, by IV. 5. a^^™^ gives the regular level surface 
of the function ZtJ'j^j (/3, h, a, t, r) fixed by the value Z^(^^_^ . 
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All the considerations concerning the dependence on the parameters t, r in the previous 
section carry over directly to a^^™^(t, r). In particular, one has 

-g^ = v^iaj^\'{t,r),t,r) , (4.20) 

where 
with 

KimM^^) ^ InFo^(m) ^ a^^'^^AC-' ({c,(n,a,r)}) > . (4.22) 

Again, we always assume that h is chosen such that dh/dt is negative. Then, from (j4.21|) . 
i;"*" > on < a < 1, with v"'' = at a = and a = 1. Also 

= Ql ^l(m) («A,fe it,r),r). (4.23) 

The derivative w.r.t. r is similarly given by ()B.3p . 
The values obey 

<5+'<a+J;^^(t,r) < l-<5+ (4.24) 

with lattice-size independent, positive and 6~^'. Again, the lower bound is automatically 
satisfied, whereas the upper bound is ensured by letting the parameter r vary, if necessary, 
in (j3.3ip (cf. Appendix B). From this it follows that the analog of (j3.32p : 

^^(i,r) > 7?+(5+) > 0, __(«+M(t,r),t) > 4(6^) > 0, (4.25) 

holds for some lattice-size independent r]f, 772". Since, furthermore, (j4.24p holds for any r 
if it already holds for r = 1, we may again set r = 1 — e, and, according to the convention 
introduced in the previous section, write C(j^^^\t) = Oi~j^^'^\t, 1 — e), etc. 

As in the last section, one may iterate this procedure of performing a decimation trans- 
formation to produce upper and lower bounds according to (j4.17p . and then fixing the value 
()4.18p of the interpolating parameter a according to (j4.19p . Assume that we choose the same 
interpolation family h at every step. Then starting from the original lattice, after n iterations 
one obtains 

Z+(/?) = 1(Za(/3) + 4-)(/3) 



2 



ni^o(m,/.,aX^^(t„,),t„)l^l/^" 

. m=l 



Z+„,({c,(n,a+^^(i„))}). (4.26) 



The discussion in subsection 13 .41 concerning the representation (j3.35p of Z\ applies equally 
well to (|4.26p . In particular, note that again the existence of the large volume limit implies 
that 

o^tr\t,r) = at^-\t,r) + 6at[r\t,r) (4.27) 
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with 6aXPit,r) as some inverse power of lattice size in the lA^"^)] oo hmit. Alter- 
natively, ()4.24p already implies that one must have a lattice-size independent contribution 
a'l;^"'\t,r) > 6+' > in (OTl) (cf. Appendix B). 

Again, either scheme (j3.33p or (j3.38p may be used to obtain (j4.26p . For the reasons 
already noted, however, the latter scheme is more convenient for our considerations. Note, 
furthermore, that the bounding coefficients c^(m) and cj'(m) in this scheme are the same for 

Za and since they do not depend on a^^f^ or a~j^^^ ■ We, therefore, adopt it in what 
follows as the common iteration scheme for Z\ and Z^: 

cm 

(4.28) 



{cf{l)} < {c,(l,ai^/,(ii))}, {c,(l,4;;Ui^))} < 
{cf{2)} < {5,(2,ag\(i2))}, {5,(2,4(^1(4))} < 



In (j4.28p and in the following, the more detailed notation h'^ and t"*" is used for the choice 
of interpolation and t-parameter values occurring in (j4.26p whenever they need be distin- 
guished from those used in the representation (j3.35p for Z\, which can, of course, be chosen 
independently. 

As indicated by the notation, even for common choice of interpolation h = and of 
all other parameters, the values of aj^^^\t,r) fixed by the requirement (j4.19p are a priori 

distinct from those of a^fl^it, r) fixed by (|3.24p . It is easily seen, however, that for sufficiently 
large lattice volume they must nearly coincide. We examine this difference more precisely 
below. 



5 The ratio Zl V^a 

We may now compare Z\ and Z\ + Z^^ ^ by means of their representations (j3.35p and (j4.26p 

on successively decimated lattices. Consider then the ratio of Za + Z^^ ^ and Za as given by 
(j4.26p and (j3.35p with common choice of interpolation h = after one decimation: 

2Z+,) (/3,/^,a+y(t+),t+) 

— ~ (1) '■^■^) 

^A(i) (^>^>«A,h(0' *) 

[ Z^,r,{[3,h,a^^\{t),t) )[ ZAa)({5,(l,<(;,^(i+))}); 
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By construction, the r.h.s. is invariant under independent variations of t and t~^. Now since 
by IV. 1 

1< 1 + ^ < 2 and 1< 1 + ^4 -r^. f < 2 , 5.3 

it follows tha1@ 

n < — z TT) < 2 . (5.4) 

Though the bounds (jS.Sp are rather crude, the resulting constraint (|5.4p is quite informative. 
First, it says that if in the equality (j3.24p . i.e. 

Za(i) {(i,h,J-l\{t),t) = Za 

one substitutes for Oi~Ahi^) wrong level surface a^'')j'*(t), the resulting discrepancy in the 
free energy per unit volume is at most 0(1/|A(^)|). Furthermore, (j5.4p constrains by how 
much a^y} {t) can differ from oi^l^ {f^) at t = t"*". From the definition (j3.19p and III.3, the 
change in ^^(1) (/^' ^' ^) under a shift 5a in a satisfies 

I 5 In (/?, h, a, t, r)| > I 5a| | A^^) | In {!) ^^^^ . (5.5) 

When combined with (|5.5p . the constraint (j5.4p . taken at general r, implies that one must 
have 

|at(;,)(t,r)-a« (t,r)|<0(-i-). (5.6) 



•|Aa)|- 



This implies that in (j3.40p . (j4.27p one has a"^\t) = a^^^^^t), i.e. any difference occurs only in 
the parts (5a^^^, Sot^^f^ that vary inversely with lattice size. Thus, in the large volume limit, 
this difference becomes unimportant if one is interested only in the computation of partition 
functions, or bulk free energies. This, however, is not the case for free energy differences such 
as the ratio (|5.ip . Indeed, any discrepancy of the size (15. 6p means that the first factor in 
(j5.2p can contribute as much as the second factor in round brackets. Thus the expression for 
the ratio of the twisted to the untwisted partition function given by (j5.2p . though exact, is 
not immediately useful for extracting this ratio on the coarser lattice. 

To address this issue one may make use of the t-parametrization invariance of (j5.2p . First 
the cancellation of the bulk energies generated in the integration from scale a to ba is made 
explicit as follows. For any given t^, choose ti in Zp^(\) (/3, a^^^^(ti), ti) so that 

Ka^l\[h),t,) = h{al^^^{tt),tt). (5.7) 

This is clearly always possible by (1332)1 and (lOSl) . and by in fact, h-f]; = 0(1/|A(i) |). 
Then (15. ip assumes the form 

V ^A ; Z^a) ({5,(l,aA,,,(ti:'^' 



(|5.4p clearly holds for general values of the parameter r, not just for the values (|3.3ip used in (|5.2p 
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We may now iterate this procedure performing (n — 1) decimation steps according to the 
scheme (I4.28p . at each step choosing tm, such that 

h{ a^;^l{U,tm) = h{ al^;^^ (t+ ), i+ ) , m = 1, . . . (n - 1) . (5.9) 
Carrying out a final n-th decimation step one obtains 

ii+?Al\ = (A/.,4y(t+),t+) ^^^^^ 

V y Z^(„)(/3,/i,ai"\(t), i) 

Z 



The argument for n = 1 (eq. ()5.2p ) above may now be apphed to (I5.11|) to conclude 

l«A??(*,0-«a(*,OI<0(^). (5.12) 

Any such discrepancy between o^}^^ {t) and Ci"f^\{t) in (jS.lip presents the same problem for 
extracting the ratio at scale as at scale ha. In this sense (jS.lip is not qualitatively different 
from the n = 1 case ()5.2p . Transferring the discrepancy to large n, however, allows a technical 
simplification as we see below. 

Next, consider (IS.lOp rewritten as 



Z+„A(^^h,a^;^\{t),t)) \ Z^(„)({c,(n,ai"i(t))} 



by use of (|4.19p . By construction (cf. (|3.24p ). a^A \(^) such that the r.h.s. in (jS.lOp . hence 

in ()5.13p . is invariant under changes in the parameter t; but note that the two a^^^-dependent 
factors in round brackets on the r.h.s. in (|5.13p are not separately invariant. If, for some 
given t, a^A^i^) larger (smaller) than a~j^^^\t), then, by IV. 5, (/?, h, a^^\{t), t) is larger 
(smaller) than .^^(„) {P,h,a^^'^\t), t) = Z'^(^„_^, and the second factor in round brackets on 
the r.h.s. of (j5.13p overestimates (underestimates) the ratio ^ /Z\. It is then natural to 

(n) 

ask whether there exist a value t = 't\ h such that 

ZU h, a^;:y;:l), tg) = . (5.14) 

Note that the graphs of a^^lii^) ^^'^ '^~Ah'^(^) ™ust intersect at t^^\- 

A unique solution to (|5.14p indeed exists as shown in Appendix C provided 

with r in (|3.3ip . An equivalent statement to (|5.15p is 
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where Aj^(^)(a,r) is defined by (|4.22p but with replaced by -^^{i)- Assume now that 

under successive decimations the coefficients c^{m) in (j4.28p evolve within the convergence 
radius of the strong coupling cluster expansion. Taking then n in (IS-lOh sufficiently large, 
we need establish inequality (j5.15| ^ only at strong coupling. Within this expansion it is a 
straightforward exercise to establish the validity of (jS.lSp , with strict inequality on any finite 
lattice. 

We summarize the above development in the following: 
V.l Consider n successive decimation steps performed according to the scheme 14.28^ . 
Assume that there is an uq such that the upper bound coefficients c^ (n) become sufficiently 
small for n > hq. 

Then the ratio of the twisted to the untwisted partition function on lattice A, of spacing 
a, has a representation on lattice A^"^, of spacing h^a and n > no, given by: 



(5.17) 



*(") _ An) / ,(n) N (r, 



where 

* (n) 

Here, the function ^-^ defined by ^3.24\ ), i-e. is the solution for a to 

^A(") t) = 2'^(„-i) , (5.19) 

and t^{^\ is defined by \5.14\) , i-e. is the solution for t to the equation 

Zl,^,{(i,Kail{t),t) = Zl,„_,, . (5.20) 



As indicated by the notation in (jS.lSp . any dependence on h must cancel in a^^"\ Indeed, 
Cauchy's form of the intermediate value theorem gives 

InZ^) ({c,(n,a)}) -In Z^^) ({ c, (n, a^^")) }) _ ^^(O ^ ^ ^^^^^ 
InZ^W ({c,(n, a)}) -In ({cj(n, «;("))}) ^aw(^)" 

for some ^ between a]^^"^ and a, and use of (j5.16p was made to obtain the last inequality. 
Setting a equal to 1 in ()5.2ip . combining with (I5.17p . and using III.3, IV. 5, gives 

> — V • 5.22 

Z^(„)({c7(n)}) 

The upper bound coefficients in (j4.28p then, which correspond to upper bounds for the 
partition functions Zp^ and \ give a lower bound for the ratio V-^aEII Setting a = 



^"For Abelian systems, comparison inequalities of the type (|5.15|l either follow from Griffith's inequalities, 
or can be approached by the same methods. All such known methods fail in the non-Abelian case. 
^^This result was first stated a long time ago in [12| . 
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in (j5.2ip . similarly yields an upper bound. Thus: 

V.2 With the same conditions as in V.l the ratio of the twisted to the untwisted partition 
function on lattice A of spacing a is hounded on lattice A^") of spacing 6"a hy: 



Z 



A(")({4(n)}) Z^(„)({c7(n)}) 



Now, the ratio of the interpolating partition functions (j4.15p and (j3.20p interpolates mono- 
tonically between the upper and lower bounds in (j5.23p since 



^ 4(i)(fe(n,a}) 
da Z^M({cj(n,a)}) 



< (5.24) 



by (|5.16p . It follows that there exist a unique value a^^"^ of a at which this ratio of the 

interpolating partition functions equals Z^^ ^ jZi^. This is a restatement of (I5.17p . but makes 
explicit the fact that this value is independent of h. In fact, it shows that all dependence 
on parametrization choices, i.e. the choice of parameters tm made in successive decimations, 
eventually cancels in a^^"^ . Indeed, the latter can depend only on the number of decimations 
n and the initial coupling /3, since this is all the upper and lower bounds in (j5.23p depend 
on. This, in retrospect, is as expected, since all bulk free-energy contributions depending on 
such choices were canceled in finally arriving at ()5.17p . but V.2 makes it manifest. 

()5.23p was obtained as a corollary of (j5.17p . An alternative approach would be to proceed 
in the reverse direction, i.e. establish (I5.23P directly, from which (I5.17|) would follow by 
interpolation between the upper and lower bounds as in the previous paragraph. In other 
words, follow also in the case of the ratio of the partition functions the approach followed 
separately for the untwisted and twisted partition functions in the previous sections. This is 
further discussed in Appendix C. 



6 Confinement 
6.1 Order parameters 

The vortex free energy ^ is defined by the ratio of partition functions considered in the 
previous section: 

exp(-F<-'(ffl) = gfi . (6.1) 

It represents the free energy cost for adding a vortex to the vacuum, the Z(l) flux of the 
inserted vortex being rendered stable by wrapping around the toroidal lattice. As has been 
discussed in the literature, all possible phases of gauge theory (Higgs, Coulomb, or confine- 
ment) can be characterized by the behavior of (16. ip as one lets the lattice become large. In 
particular, having taken the vortex to wind through the lattice in the directions k = 3, . . . , d, 
a confining phase is signaled by the asymptotic behavior 

Fi-)(/3)~Lexp(-a(/3)|A|) , (6.2) 
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where L = nK^i,2 and A = L1L2. (j6.2p represents exponential spreading of the flux 
introduced by the twist on the set V in the transverse directions (creation of mass gap), with 
(t(/?) giving the exact string tension. Note that, according to (|6.2p . — > as |A| — > 00 
in any power-law fashion, i.e. one has 'condensation' of the vortex flux. The behavior ()6.2p 
is dictated by physical reasoning [9], [10], and explicitly realized within the strong coupling 
expansion. As such free energies differences are generally notoriously difficult to measure 
accurately, demonstration of the behavior ()6.2p by numerical simulations at large /3's has 



been achieved only relatively recently [13 
The Z{2) Fourier transform of (j6.ip 



exp(-FS'(,3)) = i(l-^) (6,3) 

gives the corresponding dual (w.r.t. the gauge group center) order parameter, the color 
electric free energy. (j6.ip and (j6.3p are ideal pure long-range order parameters. They do not 
suffer from the physically irrelevant but technically quite bothersome complications, such as 
loss of translational invariance, or mass renormalization and other short range contributions, 
that arise from the explicit introduction of external sources. Such external current sources are 
introduced in the definition of the Wilson and t'Hooft loops. Furthermore, the behavior of 
the latter can be bounded by that of (16. ip and (j6.3p |15j . In particular, the following relation 
holds. Let C be a rectangular loop of minimal area S lying in a 2-dimensional [12]-plane. 
Then [l5]: 

W[C])^ < , (6.4) 



where 14^ [C] = Xi/2 \ Ylbec ^bj is the usual Wilson loop observable. It follows from (|6.4p that 
confining behavior (j6.2p of the vortex free energy implies confining behavior ('area-law') for 
the Wilson loop. 



6.2 Strong coupling cluster expansion and confinement 



ficients Cj(n,a^^\{t)) in our partition function representations (j3.35p and (j4.26p . This flow 



We now return to our considerations at the end of section [3] regarding the flow of the coef 

^A,h 

is bounded from above by that of the MK coefficients cj (m) regardless of the specific value 
assumed by the aj("/](im)'s at each decimation step (cf (I3.42P ). Furthermore, by explicit 
evaluation under the iteration rules (j2.19p - (j2.22p . one finds that c^(n) — > as n — > cx) for 
any initial (3, provided d < 4. Thus, given any initial (3, one may always take the number 
of iterations n large enough so that the coefficients {n) become small enough to be within 
the region of convergence of the strong coupling expansion. Then by V.l: 

The vortex free energy may then be evaluated in terms of the coefficients Cj{n, a^^"^) directly 
on lattice A(") of spacing b"'a within a convergent strong coupling polymer expansion. 
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Recall that, in the pure lattice gauge theory context, a polymer is a set Y of connected 
plaquettes containing no 'free' bond, i.e. no bond belonging to only one plaquette in Y (see 
e.g. [16j). The activity of a polymer Y is defined by 

^(^)= / li dUb l[9piU,n), (6.6) 

where 

gp{U,n) ='^djCj{n,a\'-''^)xj{Up) . (6.7) 
The polymer expansion is then 

lnZ^(„)= «W n (6-8) 

where the sum is over all linked clusters of polymers in A^"^ each cluster X consisting of a 
connected set of polymers Yi, i = 1, . . . ,kx with multiplicities nj. The combinatorial factor 
a{X) is given by 

a{X)= ^(-l)'(^), (6.9) 

G{X) 

where the sum is over all connected graphs on X (full set (including multiplicities) {Yi} as 
vertices with a line connecting overlapping polymers) and 1{G) is the number of lines in the 
graph. 

In the case of -^^(2), the presence of the flux enters the activities (16. 6p through the re- 
placement (j4.2p . We denote the resulting activities by z^~\Y). This replacement does not 
affect polymers that are wholly contained in a simply connected part of A^"^ since, in this 
case, the flux can be removed by a change of variables in the integrals in (j6.6|) . Only clusters 
that contain at least one non-simply connected polymer forming a topologically non-trivially 
closed surface can be affected. Thus, one has 

in4-i-inz^(„)= J2 «w ( n ^^^Ky^r- n <y^r] > (e-io) 

where the sum is only over all such topologically nontrivial linked clusters, the contribution 
of all other clusters canceling in the difference. The minimal cluster of this type consists of 
a single polymer which is a 2-dimensional plane IT : =const., /i = 3, . . . , d on A^"\ thus of 

size = L^i^L^2\ ai^d activity 
^(n) = 2^ Cj{n,al ') =cii2{n,a^>) [1+ \^ ] • (6.11) 



half — int. half — int. 

j>l/2 j>3/2 



ci/2(n, a*("))^ 



(Note that the terms from the higher representations in (16.1ip become utterly negligible 
in the large volume limit.) There are L^") = nK^i,2 -^^"^ such minimal clusters giving the 
leading contribution in (j6.10p . This leading contribution is thus seen to give the confining 
behavior (|6.2p . Nonleading contributions come from nonminimal clusters consisting of 11 with 
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or without 'decorations', and additional polymers touching 11. Such corrections have been 
evaluated in terms of the character expansion coefficients (the Cj(n, a^^"^)'s in our case) to 
quite high order [16]. They can be shown to exponentiate, so that 

exp(-6-A A) 



(6.12) 



with 



1 



n{n,at^)+0[{c^+y^/ci/2)'''^')+0{n/A^^)) 



where [T6 



K(n, a 



lnci/2(n,a^*^"^) - 4ci/2(n, a^''^* + 8ci/2(n, a]^^"^'' + 



*(n)N6 



(6.13) 



(6.14) 



By the convergence of the expansion [T7], the large volume limit exists and is given given by 
a = K(n, a*("))/&^'^, where a*^"^^ is the lattice independent part of a^^"-* (cf. ()3.40p ). 

The number of iterations n in the above expressions is taken large enough so that, given 
some initial j3 on A, the resulting c^(n) are within the expansion convergence regime, and 
one can write the representation ()6.5p by V.l. This implies the existence of a scale, a point 
to which we return below. Otherwise, n is arbitrary. By construction, our procedure is such 
that the ratio (16. ip is reproduced under successive decimations. Thus, given (16. 5p at some n, 
suppose one performs one more decimation to lattice A^""*"^). The condition that determines 
^*(n+i) g^i^]^ ih&i (|6.5p is preserved is 
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/ I 1 *(n+l) 
KA("'+l,aA 



(6.15) 



which then results in constant string tension a\ under successive decimations. Using (|6.14p 
with (j3.12p and (j2.37p . it is an easy exercise to solve (j6.15p . at least to leading approximation, 
for The i-parameter value t^i^~^^^ this a*^'^^^^ corresponds to can then also be easily 

obtained, if desired, from 

/i(a*W,i)lnFo^(n + l)|A("+i)| +lnZ^(„+i) ({ c,(n + 1, a* ("+i)) = In Z^M ({ c,(n, a* W)) , 

with lnZ^(„), lnZ^(„+i) given by (j6.8p - in fact, to leading approximation, lnZ^(„+i) can be 
ignored. Note that this amounts to replacing the set of the two equations (|5.19p - (|5.2Up in 
V.l by their ratio and one of them. This is indeed the most convenient procedure once (j6.5p 
has been achieved. 



6.3 String tension and asymptotic freedom 

a*(")) is the string tension in lattice units of lattice A^"). It is a complicated, but well- 
defined function of the original coupling /3 = ^/g^ defined on lattice A, eq. (j2.ip (cf. remarks 
following (j5.24p ). We write 

K(n,a*(")) = a(n,5). (6.16) 
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In dimensional units the asymptotic string tension in d = 4 (j6.13p - (j6.14p is then 

= \iy{9)- (6.18) 

Here, as remarked above, n is assumed greater than some required smallest n{g). This 
(dynamically generated) physical scale, or some chosen multiple of it, is the only parameter 
in the theory. Fixing it specifies how the coupling g must vary with changes of the (unphysical) 
lattice spacing a. 

It is convenient, and customary, to introduce a fixed scale Aq serving as an arbitrary unit 
of physical scales. Setting 

Ao^i = a6", (6.19) 
determines the lattice spacing a such that it takes n steps to reach length scale I/Aq: 

n = -^ln^. (6.20) 
in a Aq 

Fixing the string tension, given in units of Aq: 

a = kKl, (6.21) 

implies 

a(n,g)=k (6.22) 

for some constant k. ()6.22p specifies the dependence of the bare coupling g on n, hence, 
through (I6.20|) . the dependence on the lattice spacing a. It gives then the value g[a) specified 
by the value of the string tension. (This is, of course, equivalent to fixing (|6.18p directly.) 

Since 

a{n + l,g + Agf^ - a{n, gf^ = b [ a{n, g + Ag)'/^ - a{n, g)^'^ ] + {h- l)a{n, gf'' 
and Aa = — (6 — l)a/b for An = 1, one has from (j6.22p : 

/\^/& Ag nr 

-—(a—) = Va 6.23 
Ag Aa 

If {aAg/ Aa) = P{g), the 'beta- function', is known, (|6.23p can be integrated directly for a{g). 
(This introduces a dimensional integration constant which can serve as the scale Aq). This 
is in fact the familiar textbook argument were one assumes the existence of a string tension 
so as to get (|6.23p . in which the standard weak coupling perturbative expression for the beta 
function is then used. 

For us, however, the existence of a non-zero string tension is the outcome of the process 
of successive decimations to coarser scales as developed above. This process embodies all 
relevant information in the theory. In particular, it also supplies the specification of the 
function g{a). 

One can indeed construct the function g{a) directly as follows: 
(i) Starting with some initial value of /3 = 4/^^ perform successive decimations following the 
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flow into the strong coupling regime with resulting string tension K(n,a*^"^^), eq. (j6.14|) . at 
some n = riQ. Let k denote the value of this string tension. The corresponding value of the 
lattice spacing oq is given by (|6.20p . and g = g{aQ). 

(ii) Fix the string tension as in (16.22p . This is then satisfied at no, g{ao). 

(iii) Vary g away from (7(00) to determine g such that, under successive decimations following 
the flow into the strong coupling regime, the resulting string tension satisfies (j6.22p for n = 
no + 1. 

(iv) Repeat (iii) for n = no + 2, no + 3, • • • , no — 1, • • •• 

This provides the functional relation g{a). In particular, for b = 2, it gives the sequence of 
values (7(00/2'), / = 1, 2, . . ., starting from some value q(an)F^ 

Note that, according to (i) above, the number of decimations no at which one chooses 
to apply V.l to obtain ()6.5p . ()6.14p amounts to fixing the string tension. This is the only 
physical parameter in the theory. A specification of Ao is a specification of the value g{ao) 
at spacing ao, which is a convention of no physical import. 

One then has in principle a constructive method for obtaining g{a) by a sequence of simple 
algebraic operations. This is the coupling g{a) as defined in the physical non-perturbative 
renormalization scheme specified by keeping the string tension fixed. 

A straightforward illustration of the method is provided by setting all = 1, i.e. apply 
it to the flow according to the upper bound coefficients in (I3.33p . This yields g{a) as 
given by MK decimations. We cannot apply it explicitly to the case of interest, i.e. the flow 
following the middle column coefficients in (j3.38p . since we do not determine them explicitly 
in this paper. The qualitative features at strong and weak coupling, however, are readily 
discernible. 

At strong coupling, i.e. small initial /?, the number of decimations needed to reach a 
given string tension is of order unity, i.e. the lattice spacing a is large: a = 0{Aq^), and 
one is very far from any continuum limit. Successive decimations, by construction, reproduce 
the behavior seen within the strong coupling expansion, and the familiar strong coupling 
variation given by P{g) ~ glng is the result, as can be checked by a short computation. 

The opposite limit of large initial P corresponds to large number of decimations, hence 
a <ti Aq~^. Indeed, recall that 51 = is a fixed point of the decimations. Hence, for /3 — > 00, 
one necessarily has n ^ 00 in order for, say, the leading upper bound coefficient c^gl^) 
to reach any prescribed value < 1. Thus, aAo — > 0. Note that this limit is well-defined 
by construction since everything is bounded and continuous under successive decimations. 
Asymptotic freedom, i.e. the statement that g{a) — > as a ^ 0, is then a direct qualitative 
consequence of the flow produced by the decimations. 

It is instructive to examine the actual manner in which g{a) under the upper bound 
decimations, i.e. the Cj'(m)'s in (|3.38p . Comparing two g values that differ by one decimation 
step (5 = 2), one finds 

-Yn = -27^ + In 2 + 0(g2) (6.24) 

for sufficiently small g{a). The constant 69 = (1 — l/6^)/(241n6) underestimates the value 
ll/247r2 obtained in a continuum perturbative calculation by only about 3%. 

^^This is the analog in the present context of the 'staircase' procedure in [18) . 
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The actual flow (middle column in (j3.38p ) is faster, corresponding to somewhat larger 6o- 
According to RG lore, a beta-function defined by other means, such as fixing some renormal- 
ized coupling within weak coupling perturbation theory, should coincide, in its universal first 
two terms, with that defined by the above physical non-perturbative scheme. This, however, 
is outside the scope of, and not of direct relevance for the main argument in this paper. 

To reiterate, the above procedure completely specifies the dependence g{a) in the physical 
renormalization scheme defined by keeping the string tension fixed, and this dependence is 
necessarily such that g{a) — > as a — > 0. 

7 Concluding remarks 

In summary, we obtained a representation of the vortex free energy, originally defined on a 
lattice of spacing a, in terms of partition functions on a lattice of spacing ab"^. The effective 
action in this representation is bounded by the corresponding effective action resulting from 
potential moving decimations (MK decimations) from spacing a to spacing aft". The latter 
are explicitly computable. Confining behavior is the result, starting from any initial coupling 
g on spacing a, by taking the number of decimation n large enough. 

It is worth remarking again that in an approach based on RG decimations the fact that 
the only parameter in the theory is a physical scale emerges in a natural way. Picking a 
number of decimations can be related to fixing the string tension. That this can be done only 
after fiowing into the strong coupling regime reflects the fact that this dynamically generated 
scale is an 'IR effect'. The coupling g{a) is completely determined in its dependence on a 
once the string tension is fixed. In particular, g{a) — > as a — > 0. Note that this implies 
that there is no physically meaningful or unambiguous way of non-perturbatively viewing the 
short distance regime independently of the long distance regime. Computation of all physical 
observable quantities in the theory must then give a multiple of the string tension or a pure 
number. In the absence of other interactions, this scale provides the unit of length; there are 
in fact no free parameters El 

There is a variety of other results related to the approach in this paper that could not 
be included here. We note, in particular, that the same procedure can be immediately 
transcribed to the Heisenberg SU (2) spin model. Also, apart from analytical results, the 
considerations in this paper may be combined with Monte Carlo RG techniques to constrain 
the numerical construction of improved actions at different scales, a subject of perennial 
interest to the practicing lattice gauge theorist. We hope to report on these matters elsewhere. 

This research was partially supported by NSF grant NSF-PIIY-0555693. 

A Appendix 

In this appendix we obtain the lower and upper bounds II. 1, IV. 2, and III.l, III. 2, IV. 3, IV.4, 
and also IV. 1. 

^''This is part of the meaning of the common saying "QCD is the perfect theory". 
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§1. To prove II. 1 take the lattice A^'") to have length L^""^ = 2"^'', with integer 
in each direction = l,---,d, and torus topology (periodic boundary conditions) in all 
directions. Choose a hyperplane vrj^ without sites perpendicular to, say, the x^-direction 
and bisecting A("\ so that A^") = A^^ U A^'' with A^^ = 7^[A^"^], where TZ denotes reflection 
in vTi. Dropping now the terms in (j2.9p coming from all the plaquettes bisected by tti, all 
non-negative by reflection positivity, gives 



Z A(n) {{cj (n)}) = J d^x\(„) >(^j 



(A.l) 



Note that the 'half-lattice' A^^ has a boundary with resulting free boundary conditions for 
in the x^-direction. [d/jP^^^-^ still has periodic boundary conditions in all the other 

directions.) Next take a plane vr^ bisecting A^^ so that A^^ 
7^'[Af^'], where 7^' dem 
the r.h.s. of (|A.ip gives 



A^")'uAi;^)' with A^^^' 



TZ' [A^^ '] , where TZ' denotes reflection in vr^ . Removing now all plaquettes bisected by it' on 



Proceeding in this manner one arrives at 



> 



n) 



> 



dfA (n) 



(A.2) 



(A.3) 



In (|A.3p A^"^(l) denotes the lattice resulting from A^") by reducing its extent to one lattice 
spacing in the x^-direction, and is computed with free boundary conditions in the 

the x^-direction. One next chooses a hyperplane normal to one of the remaining directions 
bisecting A^^(l). Iterating this procedure in each successive direction, one eventually arrives 

|A 

> (/ d^tl^ , (A.4) 
where h denotes a hypercube, i.e. a subset of L^'") with 2'' sites. But, as it is easily seen. 



df,l>[l + Y.d'cj{nf 



(A.5) 



Inserting (jA.SP in ()A.4p completes the proof of II. 1. Note that, since the number of bonds is 
larger than the number of plaquettes, |A^")| in ()A.4p may be replaced by the latter. 

§2. To obtain IV.2 let 



1+n 



l+Hexp [A^~\Up,n)-A{Up,n) 



(A.6) 



'^'^ Actually a pair of hyperplanes because of the toroidal lattice topology. Following common practice, when 
a pair is actually meant, it will not be explicitly pointed out for brevity. 
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Then 



^i„)({c,(n)}) 



J dn\{n) Py Py, 



(A.7) 
(A.8) 



where in the second equahty V' is any other coclosed plaquette set homologous to V (cf. 
remarks preceding eq. (|4.2p ). As in the text, we take V to wind around the lattice in the 

{-) 

and hence -^^(n) , are clearly reflection positive in hyperplanes perpendicular to any one of the 
direction fi ^ 1,2. (jA.SP makes it clear that the measure in is also reflection positive in 



directions perpendicular to x^- and x^-directions. The measures defined by Zj^(n) and Z 



planes normal to /i = 1 or 



2: simply take V' 
■+ 



TZ[V] where TZ denotes reflection in such a 



plane. Thus, the measure in Z^^^^ possesses RP in all planes. IV. 2 (i) is a trivial consequence 



of this fact. 

To establish IV. 2 (ii) let a hyperplane vr normal to the x^-axis bisect the lattice A^'") so 
that A(") = A^"^ U A^. Take V C A^, and V = TZ[V] C A^^^ Dropping the terms in (IXSl) 
coming from all plaquettes bisected by vri, all of which are non-negative by RP, gives then 



(A.9) 



But A^ ^ has a boundary with resulting free boundary conditions for („) in the x -direction 
This implies that 



n) ; 



(A.IO) 



since one may, by a shift of integration variables, move the location of the twist-carrying set 

(n) 

V to this boundary and, by virtue of the free boundary conditions there, remove it from A^ . 
The rest of the argument then proceeds exactly as in §1 above. One thus arrives at (j4.10p . 

§3. The lower bounds III. 2 and IV. 4 are simple corollaries of II. 1 and IV. 2. Let |A(")| denote 
the number of plaquettes in lattice A^"). One has 



since |A("-i)| > |A(")|. So, by Il.l(ii) 
Z^(n-l) ({cj(n - 1)}) > 
> 



|A( 



"-1)1 



> 



|A(")| 



(A.ll) 



n 



|A(")| 



dU 



AW n [i+E*^. 

pGA(") jVO 
L 



■] Cj {n)xj{Up 



^A(")({SH}) 



(A.12) 



which gives III. 2. 
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Similarly, from (jA.lip . IV.2(ii) and using (|i3D 



^A{n-i)({c,(n-l)}) 



> 



(A.13) 



= Z+„,{{cf{n)}), 

which gives IV.4. Note that, as it clear from ([XTT]) . both (|A?T2]) and ([AlS]) are strict 
inequalities except in the trivial case where all Cj{n — 1) vanish. 

§4. To obtain III.l consider the decimation operation on the partition function Z^(m) {{cj (m)}) 
of the form (j2.9p . Let Q^^j denote the set of all [//i/]-plaquettes in A^"^). Consider the set 
of all 3-cells in A^"^) of side length b (in units of lattice spacing) in the ^-direction and unit 
side length in the fj,- and z/-directions. The basic moving operation consists of moving the 
{b — 1) interior [/ii^]-plaquettes of each cell along the positive ^-direction to the location of 
the [/i;y]-plaquette in the boundary of the cell (Figured]). Then Q[^jy] = Qj^,^] U where 

Q\pv] moved [/xi/]-plaquettes and Qy^^-^ the set of all 'receiving' [/xi/]-plaquettes 

on the 3-cell boundaries. The action of the receiving boundary plaquettes is renormalized as 
in (I2.16|) . Given Z^(m) ({cj (m)}) with corresponding action Ap{U,m), eq. (|2.24p . define the 
action 



^A(m) (U, m, ^) 



J2 ApiUp,m)+ J2 HCo-l)ApiU, 
PGA('") peg 



pi 



m 



m 



peg 



interpolating between the action before = 0) and after (^ = 1) the move. Then 

— Z^(™)(m,Ok=o = expA^(™)(?7,m,^) 1^=0 

= (Co-1) E (^^'(^^""^))o,A("^),«=o- ^ {MUp.m) 



(A.14) 



(A.15) 



0, A(™),g=o' 



where < — >q ^{m) c denotes the unweighted expectation with measure defined by (IXTIl) . 
Now, given a p E notes that (^Ap{Up, "i)^^ ^^^^ ? o (|A.15P does not vary with p 

along the K-direction i.e. along the direction of the move. Hence, if 



(jA.lSP gives 



Co 



(A.16) 
(A.17) 



since there are (6 — 1) plaquettes in Qj^j^j for each plaquette in Qf^„^i- But 



( 



^^A('") 



E (Co-l)^p(f^p,"^)- E Ap{\Jp,m) 



0,A('"),5 



> . 



(A.18) 
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()A.17P - ()A.18P then imply that, with condition ()A.16p . Z^(,„)(m, ancj^ In Zy^(m) (m, are 
increasing convex functions of ^ on < ^ < 1. Thus, ^^(m) (m, 0) < ^^(m) {m,l). 

A complete decimation A^"*) is performed by repeating the above basic moving- 

renormalization step in each of the available normal directions for each possible plaquette 
orientation as described in section \27\] One need only observe that, in carrying out each such 
successive step, translation invariance at ^ = for a given plaquette along the direction of the 
move holds regardless of any previous moves performed along other directions as schematically 
depicted in Figure [2l Each step then with Co = ^ results in a further upper bound on Z^(m) . 



Figure 2: Translational invariance along undecimated direction x'' at ^ = (see text) after 
decimation along x^. 

The completion of all the moving-renormalization steps results into ()3.4p at r = 1, i.e. 

Z^(,„,({c,(7n)}) < Fo''(m + l)l^""^'^IZA(™+i)({cf(m + l,l)}) , (A.19) 

with F^{m + l), cf{m + l, 1) given by (|239|) - ([2:22|) at C = Co"'"^^ = b^'^-'^\ Furthermore, it 
follows from Il.l(i), (jS.Sp and the fact that < c^(m + 1, 1) < 1 that replacing c^(m + 1, 1) 
by c'j {m + l,r) in the r.h.s. in ()A.19P gives a decreasing function in r on < r < 1. This 
completes the proof of III.l. 

The proof of III.l can be used essentially unaltered to obtain IV. 3. Introducing the analog 
of ()A.14p for the action (j4.4p in the presence of flux: 

A^MliU,m,0= E AirHUp,m)+ ^ ^ (^o - 1) m) - ^ ^A(-\Up,m), 

(A.20) 

^^Indeed, the second derivative of In Zj^(m) (m, ^) w.r.t. ^ is also positive by the usual convexity of the free 
energy, a statement that follows generally by an application of Holder's inequality. 
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one has 

p&Qt 1 pSQ, , 

(A.21) 

where < — >p ^ denotes the unweighted expectation with action (|A.20p . The observation 
that, having moved plaquettes in certain directions, there is stiU translational invariance at 
^ = for plaquettes to be moved in the remaining undecimated directions (Figure [2]) holds 
also in the presence of the flux. This is a consequence of the basic property (cf. section 
that ^^(ii) does not depend on the location but only the homology class of V. One may, for 
example, always bring the set V by a change of variables to occupy the exact same location 

with respect to the plaquette p in each expectation (^Ap \Up, "^)^^ ^^^^ g o (1A.21|) . Thus 

it is again the case that each such expectation in (IA.2ip does not vary with p £ Qf^,y] along 
the undecimated directions normal to p. Hence, when ()A.16P holds. 



("^' Ok=o = , whereas ^^7^ 



which gives the analog of ()A.19p in the presence of flux: 

47i)({9M}) < F,^{m + l)\^''^^''\z[~l^,,{{cf{m + l,l)}) . (A.22) 
Combining ([02|) with (lAT9j) and IV.2(i), ([33]) then gives IV.3. 

§5. To establish IV.l, take V such that all p G V are bisected by a hyperplane vr perpendicular 
to, say, x^. Let 

S^/^ = J2 c, in) XjiUp) , = 1 + E dj c, (n) Xj (Up) (A.23) 

7=half— int. i=int. 

denote the sums over half-integer and integer representations, respectively. Now 

pev pev Qcv p(zQ pev\Q 

Q|— odd 



where the sum is over all subsets of plaquettes Q in V with odd number of plaquettes \Q\ > 1. 



Inserting (|A.24p in {Zj^(„) — ^l(i)) one has 



^A(«)-^i^i=2 E [dUM^^ n [l + Y.djcf{n)xj{Up)]llSl,^' J] S^, . (A.25) 

QCV pGA(")\v j¥o peg pev\Q 



|Q|=odd 



Since Cj{n) > 0, all j, every term in the sum in (|A.25p is manifestly non-negative by RP in 
vr, which proves IV.l. 
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B Appendix 



In this Appendix we give some simple estimates concerning the variation of the level surfaces 
a"^l{t,r) and a^^\t,r) w.r.t. the parameters t,r. 

§1. The derivative of Q;^(t, r) w.r.t. r is given by 



dr 



B^(m) (a, r) 



where 



dh{a, t) 
da 

1 d 



+ A^(™)(a,r) 



(B.l) 



1 1 ^ 



The derivative of a^^'^ {t,r,r') w.r.t. r is similarly given by 



dr 



da 



+ ^A{m) (a, r) 



(B.3) 



with 



1 1 d v+ 

InFof^(m) I AM I 9r A""' 



(^{cj(m,a,r)}) < 



(B.4) 



§2. In this Appendix we use the short-hand notation = Zj^(„){{cj{m,a,r)}) and 

"^A("i) ~ -^ACm) a, ''')}); and also, given any set of plaquettes P C AM^ we define: 



d/i?, = n n fpi^p 



(B.5) 



By translational invariance 



A^im) {a,r) 



1 



In FU{m) ^ 



^ dcj{m,a,r) 1 /" n ^rr n 



A('") 



Since 



j¥0 

^AC-") > I 'if^Mm)\p Xj{Up) > 



(B.6) 



by RP in a plane bisecting p, one has 

Aj^(n.){a,r) < 



InFo^(m) 



InFo^(m) 



dcj{m, a, r, r') 



da 



jp{m,a,r) 
da 



(B.7) 
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In the same manner one obtains 

\BA(m){a,r)\ < 



1 ,2 . dcj{m,a,r) ^ 



m] 



dr 

gp{m,a,r) 
dr 



(B. 



By translational invariance for plaquettes of the same orientation in the presence of flux, 
the same bounds can similarly be shown to hold for the quantities ^^(m) and -B^^^j : 



^+(„)(a,r) < 



1 



lnFo^( 



and 



I^A(^)(«,OI < 



1 



dgp{m,a,r) 
da 

dgp{m,a,r) 



(B.9) 
(B.IO) 



^ InFo^(m) " dr 

Note that all these upper bounds are independent of the lattice size lA^™")!. 

Let p be a fixed plaquette, and c a 3-cube having p in its boundary dc and protruding 
from p in one of the d — 2 directions normal to p. Let Q denote the set of plaquettes sharing 
a bond with dc but not belonging to dc. Then 

^fJ'A('r^)\Qudc^jXj{Up) n fp'{Up',m,a,r) 



1 



Z 



A(™) 



> 



1 
1 

Zf^{m) 



p'€QUdc\p 

A('»)\Qu9c dj Xj{Up) n fp'iUp',m, a, r) 

p'&dc\p 



df^Aim)\Qiidc djCj{m, a, r)^ 



> 



d'j Cj{m, a, r)^ 



\\fp{m,a,r)\\\Q\[l + Ei^o df Q(m,a,r)6] 



(B.ll) 



RP in each of the two planes bisecting the plaquette p shows that the terms involving plaque- 
ttes in the set Q in the first line are all positive which results in the inequality in the second 
line. Using ([B?TT]) in ([R6]) gives the lower bound 



Aj^(m){a,r) > 



d'^Cj{m, a, r)^ 



1 dcj{a,r) 
\\fp\\\Q\ \nF^{m) f^^ ~d^[TTE^^^WWw¥] 



1 



Similarly, one obtains a lower bound on \B^(m) 



1 + ^ Cj{m, a, r)^ 



(B.12) 



|S^(m)(a,r)| > 



1 



E 



dcj{a,r) (fjCj{m,a,rf 



dr '[1 + E.^o^? 



Ci[m,a,r) 



InFo^(m) 6^ dr 



1 + ^ d^ Cj{m, a, r)^ 



(B.13) 



41 



Again, these lower bounds are manifestly lattice-size independent. 

Since, as it is easily seen (apply (|A.8p ). (jB.lip holds also when ^/^^(m) is replaced by 
d/z^j"^), the r.h.s. of ()B.12|) and ()B.13P also give lower bounds for A^(„i) and |-B^(„)|, respec- 
tively. 

§3. The lower bound in (j3.30p . i.e. 



(B.14) 



is a consequence of II. 1 and (j3.24p . From the first inequality in ()A.12p (excluding the trivial 
case of all Cj vanishing) one has 



(B.15) 



whereas from (I3.24p . with the short-hand notation a)!^''il{t,r)) = a^^l, 



Combining (jB.lSp and ()B.16p gives 



1 + Y,d]cj{m,a^^l,r) 



(B.16) 



(B.17) 



which, together with (|3.17p and (13. Sp or (|3.9p . shows that it cannot be that a"J^l^{t, r) — > in 
any fashion with increasing lattice size |A("^)|. 

An explicit lower bound on is easily obtained from ()B.17p by taking (j3.12p with, for 
example, (|3.9p and ht{a) = a^^^\ s{t) > 1. Using the elementary inequality 

x"^ -l<q{x-l) , x>Q, 0<g<l, 

(jBT7|) gives 



(B.18) 



1 > aH(t,r) > T.,^od]c,{m-l) 



( (m) - 1 ) + (m) E,yo (m, r 



'2 JJl 



5' > 0. 



(B.19) 



Note that, from (iD^ll below, {F^ {m) - 1) > T.j^od'] {m - if > T.j^od]cj{m - if. 
Similar expressions can be obtained for other choices of /ij. 

The lower bound in (j4.24p is similarly seen to hold by combining IV. 2 and (j4.19p . 

To satisfy the upper bound requirement in (j3.30p . i.e. 



atl{t,r)<l-5, 



(B.20) 



it suffices to let the decimation parameter r vary, if necessary, away from unity in the domain 
()3.3ip . To see this, suppose that, performing the m-th decimation step, one finds that 



42 



= 1 - '5a™\ where, say, 6^^^ < 0(1/| A("^)). Using the bounds ^3), ([B13]) 
on Ay^(m), B^(,n), and the boundedness of the derivatives of h{a,t), in eq. (fRTI) then ffives 

da'£l{t,r)/dr>e>0 (B.21) 

for some constant 9 independent of lA^^^I for all < r < 1. Hence, for t > tm, and some ^ 
between and 1, 



r=l-^e 



dr 

< 1 - (5(e) (B.22) 

with 5(e) = e6'/2, and (|B?20]) is satisfied. 

Given (j3.30p . the bounds (jB.7p . ()B.12p and the properties of the interpolation h (cf. 
([3T71) . msTTE^i ). it follows from ([3:25]) and ([3:29]) that 

^(t,r) > r?i(5) > 0, __(ag(t,r),t) > r?2(5) > 0, (B.23) 

where jyi, 7/2 are lattice-size independent, and r in the domain (I3.3ip . 

Since the same upper and lower bounds apply to , |, the same considerations 

again show that eqs. ([4.24p and ([4.25P are always ensured to hold by letting r vary, if 
necessary, in (I3.3ip . 



C Appendix 

§1. Under the conditions in V.l, given a^^\{t) satisfying (I3.24p . one seeks a solution t = t^^\ 
to ([OO]) . i.e. 

Z+„,(/3,/i,ag(t),t) = Z+„_„ . (C.l) 

First note that, if a solution exists, it is unique since, as it is easily checked, (/3, /i, a^]j(i), t) 
is monotonic in t. 

To show that a solution exists, we proceed as follows. By ([4.19p . given any suitable 
interpolation /i, there is a function oi^^\t) such that 

(/?,/!,<?(*), t) = ^+(„-i) (C.2) 

for all allowed values of the parameter t. We simplify notation in the following by omitting 
the fixed labels n. A, and write 
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Assume that there is a t = tj at which ah{tj) > a'l^{ti). Then let to be such that 



h{al{ti), ti) = h{ahito), to) , G [ta^h] ■ 



(C.3) 



It is always possible to have such a to by virtue of (I3.32P and (I4.25|) . which imply (Figure 
[3l^a)) that to > tj and 



(^h{t) > Oft (to) > Oft (ti 



t > to 



(C.4) 



ah(t) 



tj Iq 




h{ah{t),t) 
h(at{t),t) 




h{a+{t),t) 

h{ah{t),t) 
^ t 



(a) 



to tj 

ib) 



Figure 3: Fixing to in the neighborhood of some initial t/. Distance between curves in the 
figure is greatly exagerated. 



With the notation 
we now define 

^{X,t) ^ hia,{t),t) + (1 - A) $+„)(a+(t,)) + A $+„)(«,, (t)) - , (C.6) 
and consider the equation 

^(A,t) = 0. (C.7) 

At A = eq. (|C.7p is solved by setting t = to since there, by ()C.3p . it reduces to (|C.2p 
evaluated at t/. At A = 1, (\C.7\i becomes the equation to be solved (jC.ip . A solution to (|C.7p 
then determines implicitly a function t(A) with the property t(0) = to- If this function can 
be extended on < A < 1, it gives the branch of solutions of ()C.7p through (0,to), and t(l) 
will be the solution to the original problem (|C.ip (method of imbedding or continuity [19j). 
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By the implicit function theorem, if grad ^ is continuous and {d^ /dt){0,to) ^ 0, there 
exists a branch t{X) through (0, to) on a sufficiently small interval around A = 0. One then 
extends t{X) by a standard argument. Denoting partial derivatives by subscripts following a 
comma, one has 

For sufficiently small AA then, and by the mean value theorem, one can write 

t{0 + AA) = t{0) + t,A(^AA) AA = t{0) + t^xiO) AA + 0{A\^) 

for some < ^ < f . One can then use (jC.SP to find a(0 + AA), and repeat the procedure to 
obtain t(0 + 2AA), and so on. If grad ^ is well-behaved throughout the relevant X — t domain, 
this procedure constructs the desired branch away from the initial point as long as ^'^t 7^ 
along the branch. The existence of a solution t{l) is therefore guaranteed by basic existence 
theorems (see e.g. jl9j ) if this condition is satisfied throughout the interval < A < 1. 

Now 



*,i(A,t) 



h{ah,t)^a + A A+(afe) 



+ 1 



h{oihit),t 



< 



by (jS.lSp and since h{ah,t)^t < 0. On the other hand, 
^,A(A,t) = [«^+(„,(a,(t)) - f+„)(a+(t7)) 



> 



< A < 1 



< A < 1 . 



(C.9) 



(C.IO) 



by (IC.4p and IV. 5. Thus, from (IC.Sp . t(X) is an increasing function of A, and extends to the 
solution t(l) > to- 

Conversely, if there is a tj such that ah{ti) < Q;jJ'(t/), one can find to such that ()C.3P is 
satisfied (Figure El^b)), where now to < t/ and 



(*/) > «^(*o) > ah{t) , 
Now ()C.9P remains unchanged, but 

^, a(A, t) = [ <!>+„, (a,(t)) - (a+(t,)) 



t<to. 



< 



< A < 1 



(C.ll) 



(C.12) 



by (|C.lip and IV. 5. It follows that t(A) is now a decreasing function of A, i.e. if ah,(t/) < 
a'l^{ti) at the starting point one moves backwards in t to hit the point t(l) where ah{t{l)) = 
a^(t(l)) and (jC.ip is satisfied. This concludes the demonstration of the existence of a solution 
to dSIll). 

§2. Going back to (jS.lOp . assume there is a tj such that a^^/j(t/) < a^^^\tj). Then, setting 



t = t^ 



tj, one has 



1 + 



.(-) 



Za 



Z^(n) {(3,h,af\{tj), tj) 
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.{-) 

> I 1+ ^ ' — m (c.is) 



.{-) 



A(") ({cfH} 



by IV. 5 and (j5.24p . On the other hand an upper bound is always obtained by using (j4.5p . 
Taking the coefficients Cj{n,a) in the form ()3.12p with ()3.9p . which interpolate between this 
upper bound at a = and the lower bound ()C.14p at a = 1, one obtains 

^ ^ fUMM > £ > fMiMl) (C15) 

Z„„., {{df{n)}) Za Z^,..,({cf{n)}) 

This is (j5.23p again (for the case Cj{n) = 0). It follows from ()C.15P that there exist a value 

a*j^^^^ such that ()5.17p holds. This value is unique by monotonicity from (I5.24p . 

This is an alternative way of treating the ah{ti) < case in §1 above. 

()C.14p - (jC.lSP partially implement the alternative approach to (j5.17p outlined in the last 
paragraph of section [5l One way to complete it would be to show that, for some interpolation 
h, there is at least one value tj such that a^\\{ti) < a~j^^^\tj). 



D Appendix 



From (I2.22P with integer C > 1 one has 



Fi{n + 1) = 6o,i + (1 - <5o,.) C Ci(n) + ^i(^) ' 



(D.l) 



where 



1 



* {js\ 0<s<k} h,-,lk-i 

(D.2) 

with A{j,k,l) = 1 if j,k,l form the 'angular momentum addition triangle' relation, i.e. 
I = \ j — k\, ■ ■ ■ , j + k, and otherwise. Then 

Ydiliik) = Y dj^Cj^{n) ■ ■ ■ dj^Cj^{n) ^ A(ji, j2, ^2)A(/2, ja, ^s) • • • A(4-i, Jfc, i) 

< J2 dj, Cj, (n) • • • dj, Cj, (n) dj^ dj, • • • dj^ 

{js\ 0<s<k} 

< Mn)\\'. (D.3) 



Hence, with b > 2, 



Y.dlHn + 1 



T.dl 



62 
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< 



< 



ij^O k=2 \ / ij^O 



62 



ClbHII + E ullbNII' 



k=2 



62 



[l + lb(n)||]«-l 



62 



< 



62 



CMn)\\ 1 + Mn) 



(C-l)62 



Also, from (|m]) . (ICT 

Fo(n + l)>l + ^^Io(2) = 1 + ^^E4^.H' 

> 1, 

whereas also from (|2.22p 



Fr,(n + 1) < 



l + \\g{n) 



(D.4) 

(D.5) 
(D.6) 

(D.7) 



Combining (|D.4p and (|D.6p and taking r = 1 gives (j2.3ip . 

For r 7^ 1, (|D.4p and, hence, (|2.3ip hold with the replacement 6^ — > 6^r, provided > 2. 
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